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Abstract: 

We propose a new subtraction scheme for next-to-leading order QCD calculations. Our 
scheme is based on the momentum mapping and on the splitting functions derived in 
the context of an improved parton shower formulation. Compared to standard schemes, 
the new scheme features a significantly smaller number of subtraction terms and fa- 
cilitates the matching of NLO calculations with parton showers including quantum 
interference. We provide formulae for the momentum mapping and the subtraction 
terms, and present a detailed comparison with the Catani-Seymour dipole subtraction 
for a variety of 2 — )■ 2 scattering processes. 



1 Introduction 



The Standard Model (SM) of particle physics has been tested and confirmed by a large 
variety of experimental data [1]. With the start of data taking at the Large Hadron 
Collider (LHC) in 2009 [2], the exploration of the properties and predictions of the SM 
is being extended to center-of-mass energies of several TeV. To further consolidate the 
Standard Model is an important effort per se; however, a detailed knowledge of SM- 
induced backgrounds is equally important in searches for physics beyond the Standard 
Model. For both tasks, leading-order (LO) QCD calculations of SM processes are 
not sufficient: the theoretical uncertainty at LO is substantial, and next-to-leading 
order (NLO) corrections to LHC processes are in general large and have to be included 
to match experimental accuracies (for a recent review on higher-order corrections for 
hadron collider processes see [5]). More specifically, the analysis and interpretation 
of experimental signatures at the LHC require theoretical predictions for differential 
distributions or cross sections with cuts on kinematic variables. NLO calculations 
should thus be set up in form of fully differential parton-level Monte Carlo programs. 
In recent years, a lot of effort has gone into semi- or fully-automated inclusion of 
higher-order corrections in such tools for processes at hadron colliders. Most of these, 
as e.g. MFCM [1], NLOJet++ 0, VBFNLO [B], Rocket [7], Blackhat/Sherpa [HIE], 
GOLEM/Whizard [TO], GOLEM/Madgraph [lO], and Helac [HI [El US], are now able to 
provide NLO QCD predictions for several processes at parton level. Finally, a leading 
logarithmic description of multiple soft and collinear parton emissions in the form of 
parton showers should be combined with the NLO partonic predictions; this has been 
realised in a variety of different approaches [HIISllISlinillSlIIHlEniEIlES], and 
mature program packages like MC@NLO [ISl [23l EH ES] and POWHEG [26l EZl [28] 
have been developed for LHC physics. 

An important calculational tool for the implementation of NLO QCD corrections 
in Monte Carlo style programs are subtraction schemes [29l |30l SH ESI ESI El] , which 
facilitate the treatment of infrared and collinear divergences originating from different 
phase-space contributions. Subtraction schemes introduce local counterterms which 
mimic the behaviour of the real-emission matrix element in the singular limits. The 
difference between the local counterterms and the real-emission matrix element is thus 
finite and can be integrated numerically. The infrared and collinear singularities are 
isolated by integrating the subtraction terms over the singular regions within a given 
regularisation scheme, e.g. dimensional regularisation. The soft and collinear singular- 
ities then cancel when the integrated subtraction terms are added to the virtual cross 
section. After standard UV-renormalisation, all contributions to the NLO cross section 
are finite, and the further phase-space integrations can be performed numerically by 
means of Monte Carlo techniques. 

Various general NLO subtraction schemes have been proposed (see [30 | [3T | l32 | IM]). 
and several (semi) automated implementations are available [351 EHl EH EH Ell [28] . The 
schemes differ in the phase-space mapping which relates real-emission and leading-order 
kinematics and which is needed to define the subtraction terms, and in the finite parts 
of the subtraction. Unfortunately, the schemes developed so far suffer from a rapidly 
rising number of momentum mappings needed to evaluate the subtraction terms, which 
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basically scales like for a leading order 2 ^ N processJl] With an increasing number 
of final-state particles this scaling leads to a rapidly rising number of momentum map- 
pings and subsequent re-evaluations of matrix elements. In this paper, we therefore 
propose a subtraction scheme where the number of mappings scales like A^^, thereby 
reducing the number of matrix element evaluations by a factor A^. This scaling is 
achieved by the use of a mapping prescription which, for emissions from final-state 
particles, redistributes the momenta to all non-emitting final state particles simultane- 
ously, thereby leading to one unique momentum mapping per emitter. Furthermore, the 
subtraction terms in this scheme are derived from splitting functions which have been 
proposed in the context of a parton shower including interference effects [iQl HU H2] . 
When NLO parton-level calculations and parton showers are combined, specific coun- 
terterms have to be added in order to avoid double counting of contributions which 
are included in both the NLO calculation and the parton shower The use of 

the shower splitting functions as subtraction terms reduces the number of these coun- 
terterms [19] and therefore facilitates the combination of shower and NLO calculations. 
First results obtained with the new subtraction scheme have been published in jUj. 

This paper is organised as follows. After a brief revision of the generic features of 
NLO subtraction schemes in Section 2, we shall present the setup of our new scheme in 
Section 3. We discuss the momentum mapping, the factorisation of the matrix element 
in the soft and coUinear limits, and present the subtraction terms and their integrated 
counterparts for scattering processes with at most two particles in the final state. 
In Section 4, we apply our scheme to well-known collider processes at NLO: single- 
W production, dijet production at lepton colliders, gluon-induced Higgs production 
and Higgs decay into gluons, and deep inelastic scattering. We show that, for these 
processes, the results obtained using our subtraction prescription agree with the results 
obtained within the commonly used Catani-Seymour dipole subtraction scheme [31] . 
We conclude in Section 5. Some additional useful formulae are listed in the Appendix. 



2 General structure of NLO cross sections and sub- 
traction schemes 

In this section, we shall briefly review the general structure of subtraction schemes 
for NLO cross-section calculations at colhders and set up our notation, following 
closely the notation established in the context of the Catani-Seymour dipole subtrac- 
tion scheme [31] . 

We consider a generic jet cross-section a with 

= f da"" + f da^ + f da"^, (1) 

J m J iri J m+1 

^We note that recently a constant scaling behaviour for iV-gluon final states has been achieved 
within the Madgraph framework 38i; however, this scaling behaviour relies on the symmetrisation of 
the matrix element and is therefore applicable independent of the subtraction scheme and its mapping 
prescription. Our argument concerns the general combinatoric number of mappings and the resulting 
scaling behaviour if no further process-specific simplifications are applied. 
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Figure 1: Soft /collinear factorisation: when the partons ^ and j become colhnear 
and/or parton j becomes soft, the (j)-parton matrix element factorises into a sum over 
m-parton matrix elements times a singular factor vi. 



where cr^, o"^, and denote the LO, virtual and real-emission contributions, respec- 
tively. There are m partons in the final state for the LO and virtual cross sections, 
and m + 1 partons for the real-emission contribution. After UV-renormalisation, the 
virtual and real-emission cross sections each contain infrared and collinear singulari- 
ties, which we regularise using dimensional regularisation, i.e. we work in c? = A — 2e 
dimensions so that the singularities appear as l/e^ (soft and collinear) and 1/e (soft or 
collinear) poles. These poles cancel in the sum of virtual and real contributions, but 
the individual pieces are divergent and can thus not be integrated numerically in four 
dimensions. 

In subtraction schemes one constructs local counterterms which match the be- 
haviour of the real-emission matrix element in each soft and collinear region. Sub- 
tracting these counterterms from the real-emission matrix elements and adding back 
the corresponding one-particle integrated counterparts to the virtual contribution re- 
sults in finite integrands for both the virtual correction (m-particle phase space) and 
the real contribution (m + 1-particle phase space): 



a 



NLO 



Jm J m+1 J m+1 



finite 



finite 



I da^ + f da^ + I [daf^Q - da^^^] . 

J m \_ Jl J e=0 J m+1 



(2) 



The construction of the local counterterms, collectively denoted by da^ in Eq. ([2]), 
relies on the factorisation of the real-emission matrix element in the singular (i.e. soft 
and collinear) limits (Fig. [I]) gSl SHI HZ] : 



Mra+l{{P] 



m+lj 



^ Vi{{p}m+l) ® Mm{{p}m) , 



(3) 



where M.rn+i and Aim denote (m + 1)- and m-parton matrix elements, respectively, the 
Vi are generalised splitting functions containing the singularity structure of the m + 1 
matrix element, and the symbol ® represents properly defined phase-space, spin and 
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colour convolutions. (See Section I3l2] for a more detailed discussion of the (m + 1))- 
parton matrix element factorisation.) The momenta in the (m+1)- and m-parton phase 
spaces are denoted by {p}m+i and {p}m, respectively. As M.m+i and M.^ are defined 
in terms of (m + 1)- and m-parton phase spaces, a mapping — ?■ {p}m needs 
to be introduced. The mapping must be such that four-momentum conservation as 
well as the on-shell condition for all external particles are satisfied for both the {p}rn+i 
and the momentum configurations. Different subtraction schemes differ in the 

definition of the generalised splitting functions Vi and in the mapping from (m + 1)- 
to m-parton phase space. The momentum mapping and the construction of the vi for 



our subtraction scheme will be discussed in detail in Sections 13. H 13.21 and 13.31 

Squaring the generalised splitting functions Vi and summing over all different sin- 
gular parton splittings + j yields the subtraction terms, symbolically written 
as 



da' 



Vf ® da 



B 



(4) 



where oc da^ denotes the LO cross section, and the symbol ® represents 

properly defined phase-space, spin and colour convolutions as above. 

Integrating the subtraction term da^ over the one-parton unresolved phase space, 
d^p, yields an infrared- and coUinear-singular contribution 



d^pVi 



da^ = ^Vt® da^ 



(5) 



which needs to be combined with the virtual cross section to yield a finite NLO cross 
section 



a 



NLO 



da^ 



da^ 



J m+l '- n 



da 



B 



(6) 



In this form, the NLO cross section can be integrated numerically over phase space 
using Monte Carlo methods. 

The generalisation of Eq. ([6]) to hadron collisions is straightforward and requires 
the inclusion of a further counterterm to absorb initial-state collinear singularities 
into a re-definition of the parton-distribution functions. Finally, we emphasise that 
the jet cross-section a has to be defined in a infrared-safe way by the inclusion of a 



(m+l) 



— > p'f^'' in the collinear and infrared limits. 



jet- function Fj, which satisfies F 
Both, the factorisation of initial-state collinear singularities and the inclusion of a jet- 
function are standard and are included in the more detailed version of our final formulae 
presented in Section 



3 Alternative subtraction scheme: setup 

In the scheme proposed in this paper, the NLO subtraction terms are derived from 
the splitting functions introduced in the formulation of a parton shower with quantum 
interference [101 El 112], and the momentum mapping {pj^+i — {'p}m used to define 
the subtraction terms corresponds to the inverse of the mapping introduced in the 
definition of the parton shower. 
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In Section 13.11 we shall first describe the momentum mappings, following closely 
Ref . 110] • The essential feature of our scheme is that we use a global mapping in which 
all of the partons participate, c.f. Ref. |1H] . The factorisation of the real-emission matrix 
element in the singular limits and the construction of the subtraction terms is then 
discussed in Section 13.21 The subtraction terms needed to address the NLO collider 
processes considered in this paper and their integrated counterparts are presented in 
Section 13.31 Section 13.41 finally collects all necessary formulae. While our scheme can 
be formulated in a completely general way, here we restrict ourselves to the discussion 
of processes with massless QCD partons. 

3.1 Momentum mapping 

As before, the (m + l)-parton phase-space four-vectors are denoted by pi,P2, ■■■ and m- 
parton phase-space four- vectors by pi, p2, .... Indices a, b are used to label initial-state 
particles. 

3.1.1 Final-state emissions 

Let us first treat the splitting of final-state partons. We shall briefly review the mo- 
mentum mapping used in the formulation of the parton shower in Ref. [30] and then 
define the inverse transformation needed for the subtraction scheme. 

Mapping in the parton shower 

We begin with a m-parton configuration with momenta {p}m, where one of the par- 
tons with label £ G {1, . . - m} splits. After the splitting we have a (m + l)-parton 
configuration with momenta {p}m+i- We label the daughter partons by indices i and 
j, i.e. we consider a splitting pi ^ pe + PjH For an exactly collinear splitting or the 
emission of a soft gluon with momentum pj = one has pe = pe +Pj, but away from the 
collinear /soft limit p£ ^ pe +Pj in general. Thus, in order to satisfy four-momentum 
conservation and to have all external partons before and after the splitting on their 
mass-shell, we need to take some momentum from the spectator partons which are not 
involved in the splitting. In contrast to the Catani-Seymour scheme |3T], where the 
momenta of most of the partons are left unchanged, we use a mapping where all final 
state partons participate. 

Let us now first provide the formulae which specify the momentum mapping used 
in the parton shower of Ref. We choose to leave the momenta of the initial -state 
partons unchanged. 

Pa = Pa, Pb = Pb, (7) 

and introduce the variable Q to denote the total momentum of the final-state partons 
in m-parton phase space: 

m 

Q = ^Pn = Pa+Pb- (8) 

n=l 

^Ref. [10] uses instead of the label "j" the label "m + 1" for one of the daughter partons. 



5 



Since the the momenta of the incoming partons remain the same, Q = pa + Pb- We 
now define 

= 7^- (9) 

Note that > 1. The total momenta of the two daughter partons p^ and pj are then 
parametrised as a hnear combination of and Q according to 

= Pi + P, = Xpi + ^—7^^Q. (10) 

za£ 

The two parameters A, y in Eq. f|TO|) can be be determined from energy-momentum 
conservation as ^ 

A = J{l + yf-Aaey and y = — ^ . (11) 
The parameter y is a. measure for the virtuahty of the sphtting with a maximum value 

Z/max = {y/oe- - l) (12) 

corresponding to A = 0. 

In our scheme, the mapping includes all final-state particles, and their momenta 
before and after the splitting are related by a Lorentz transformation 

= A(i:,ir)^K, n^{e,j}. (13) 

Here K is the total momentum of the final-state spectators before the splitting 

K = Q-pi, (14) 

and K is the total momentum of the final-state spectators after the splitting 

K = Q-P,. (15) 

Since each final-state spectator is changed by a Lorentz transformation, we have 

k^' = AiK^KY^K" (16) 

with the Lorentz transformation 



^ ' ' " (k + Kf ' 

In this paper we shall focus on processes with only up to two massless partons in the 
final state; then, = 1 and = 0. In this case, an alternative representation of 
the Lorentz transformation has to be introduced which is well defined when = 0: 

\K-n J yK-n 

where n and n are light-like vectors in the Q-pe plane with n-n = 1 and [p^-n/pi-n) < 
(Q-n/Q-n). 
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Let us briefly comment on the flavour structure of the splitting. The flavours 
/ G {g, u,u,d,d, . . .} of the spectator partons remain unchanged 

fn = fn, n^{i,j}, (19) 

while the flavour of the mother parton fi obeys 

fi + fj = fi, (20) 

where we use the notation of adding flavours as in q + g = q, q + q = g etc. Thus if 
the mother parton £ is a quark/ ant iquark, then {fi/fj) = {q/q,g), and if the mother 
parton £ is a gluon, then {fi,fj) = {g,g), which corresponds to g ^ g g splitting, or 
any choice of quark/ ant iquark flavours {fe,fj) = {q,q), which corresponds to g qq 
splitting. 



Mapping in the subtraction scheme 

There is an inverse of the transformation Eq. f|T3|) . which maps the (m + l)-parton 
momenta to the m-parton momenta. It is this inverse transformation which is needed 
to determine the subtraction terms in a NLO calculation. We thus start with a {p}m+i 
configuration and determine the {p}m configuration by combining two final state par- 
tons Pi + pj -> pi. 

One can determine pe by rearranging Eq. ( Iiup 

Pe = Upi+P,)- ^~^^^ Q, (21) 

where we now have to express the momentum Q and the parameters a^, A, y defined in 
Eqs. ([8]), ([9]) and (|TT|) . respectively, in terms of the {p}m+i- From Eqs. ([7]) and ([8]) we 
have 

m+l 

Q = Pa+Pb = ^ Pn, (22) 
n=l 

while y and a£ can alternatively be written as 

y = K and ap = ^ , (23) 

with Pi = Pe + Pj (|TOi) . With y, ae as given in Eq. (l23l) . the parameter A follows from 
Eq. (dH). 

From K = Q — Pi and K = Q—P^ we can finally obtain the Lorentz transformation 
which takes the spectator partons from the (m+l)-parton to the m-parton phase space: 

p^^ = A{K,kr,pl n^{i,a,b}, (24) 

where A{K, K)^^, is given by Eq. f lT7|) with K and K interchanged. For = 1, 
i.e. {p}m = {pi,Pm}, the mapping is particularly simple and reads 

Pi = iPe + Pj-yQ), Pm = • (25) 

1 - 1/ 1 - y 
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The transformation of the flavours is similar to the case of parton splitting. The 
flavour of the mother parton fg is given by 

fi = fi + fj, (26) 

with the rule of adding flavours, q + g = q and q + q = g, while the flavours of the 
spectators remain unchanged 

fn = L n^{i,j}. (27) 

3.1.2 Initial-state emission 

Let us now turn to the description of the splitting of an initial-state parton into an 
initial- and final-state parton, and the combination of an initial-state and a final- 
state parton into an initial-state parton. As before we first review the momentum 
mapping defined for the parton shower evolution in Ref. |10] and then present the 
inverse mapping needed for the NLO subtraction scheme. 

Mapping in the parton shower 

We take the initial state partons to be on-shell with zero transverse momentum, i.e. 

pI=pI=pI= Pi = (28) 

and, in general. 

Pa = VaPA, Pb = VbPB] 

Pa = VaPA, Pb = VbPB- (29) 

Here pa and pb are the momenta of the incoming hadrons, which are taken massless 
p'a = p'b = and 2pA ■ Pb = s, where s denotes the hadronic center-of-mass energy. 
The 7]a, rjf, are the usual momentum fractions of the hadrons A and B carried by the 
partons a and b, respectively. 

We now consider the splitting of an initial-state parton, say parton a, into a new 
initial-state parton and a final-state parton, pa — )■ pa + Pj- Here, the splitting is to be 
understood in the sense of backward evolution, i.e. the evolution going forward in time 
is pa ^ Pa + Pj, where parton a enters the hard interaction. Given {p}m and pj we 
need to define how to obtain {p}m+i- As in the case of final-state splitting, we cannot 
have on-shell partons with pa = Pa + Pj away from the soft /collinear limit, so we need 
a momentum mapping which redistributes some momenta from the spectator partons. 

To define the momentum mapping we first choose to keep the momentum fraction 
of parton b unchanged 

Vb = Vb- (30) 

The momentum fraction fja will be determined from the momentum pj. To proceed, 
first consider the momenta of the final-state spectator partons after the splitting. They 
are related to the momenta before the splitting by a Lorentz transformation 

p>;, = A{k,Kr,pl nG{l,---,m} and n^j. (31) 
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Here A(i^, is the transformation specified in Eq. (fT7|) . and K and K are the total 
momenta of the final state spectators before and after the splitting, 

K = Pa+Pb 

K = Pa+Pb-Pj = Q -Pj, Q= Pa+Pb- (32) 

To determine fja one uses the fact that K and K are related through the Lorentz 
transformation Eq. ([IT]), k^' = A{k,KY^K'', so that K'^ = K'^. This condition 
determines as 

^^ = !hJM±lJ]hPi:A<i. (33) 

ribS- 2pA-Pj 

As discussed for the final-state splitting in Sec. 13.1.1^ the flavours of the spectator 
partons remain unchanged 

fn = fn, {a,j}, (34) 

while the flavours of the daughter partons fa and fj obey 

fa + fj = fa. (35) 

Mapping in the subtraction scheme 

To obtain the inverse mapping which takes us from {p}m+i to {p}m and which specifies 
the NLO subtraction terms, we first determine r]a from the condition K"^ = K"^ as 

_ fla7]bS-2 {fjaPA ■ Pj + r]bPB ' Pj) .^rX 

VbS 



which leads to 



Pa={l-^\Pa. (37) 
Pa-Pb 



As before, we choose to set rj^ = fj^. The final-state spectator partons {p}m are deter- 
mined through the inverse of Eq. (13T]) : 

p^^ = A{K,kr,pl nG{l,---,m} and n^j, (38) 

where A{K, K^jj is given by Eq. fll7p with K ^ K. The total momentum of the final- 
state partons in the m-particle configuration can be calculated from K = pa + Pb = 
VaPA + VbPB using Eq. ([36]). 

The flavours of the partons obey the usual relations 

fn = fn, {a,j} 

fa = fa + fj. (39) 

3.2 Factorisation in the soft and collinear limits 

We now proceed to derive the subtractions terms in Eq. ([!]) from the factorisation 
of the matrix-element in the soft and collinear limits, c.f. Eq. ([3]). If partons I and j 
are collinear, we have 

\M{{p, /Wi)) = \M,{{p, /}™+i)) , (40) 
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where the partial amphtude \J^e{{p, f}m+i)) can be expressed in terms of a m-parton 
amphtude times a sphtting function as specified below. Here we follow the notation of 
Refs. [211 Sn], where a QCD amplitude is written as a vector |A^) in colour and spin 
space. In the case that parton j becomes soft, the full amplitude is given by a sum of 
the partial amplitudes, 

\M{{p, = J2 /Wi)) • (41) 

e 

In the soft and coUinear limits, the partial amplitudes \J^e{{p, f}m+i)) take the fac- 
torised form 

\Me{{p, fWi)) = tlife ^fe + /,) V,\{p, \M{{p, /}„)) , (42) 

where V^{{p, f}m+i) is an operator acting on the spin part of the amplitude, while the 
operator tl{fe fe + fj) acts on the colour part. 

The m-parton amplitude |A^({p, f}m)) is evaluated at momenta and flavours {p, f}m 
determined from {p, /}m+i according to the transformations specified in Sections 13.1.11 
and 13.1.21 for final- and initial-state emitters, respectively. The matrix elements of the 
spin-dependent splitting operator can be expressed as 

{{s}m+l\V^H{p, f}m+l)\{s}rn) = JJ 5s„,s,A Vi{{pJ}m+l, Sj , Se, Se) . (43) 

The splitting amplitudes ve in Eq. (H3l) have been derived in Ref. [IQ] from QCD vertices. 
In general, they depend on the type of partons involved in the splitting. The results 
are collected in Appendix \M In the special case that parton j is a gluon and that pj 
is soft, or soft and collinear with pi, the sphtting amplitudes simphfy and are given by 
the eikonal approximation, 

vf"{{p, fU+u Si, s,) = v^sSs„s, <P^^'^^Qy-P' , (44) 

Pj ■ Pi 

where Q = Q is the total momentum of the final state partons, see Eq. ([8]). 

To construct the NLO subtraction terms, we consider the square of the amplitude 
in the singular limits, (A^({p, /}m+i)|A^({p, /}m,+i))- There are two kinds of contri- 
butions: the direct terms which correspond to the amplitude for a parton £ to split 
times the complex conjugate amplitude for the same parton to split, and the interfer- 
ence terms where parton j is emitted from parton C. in the amplitude and from parton 
k ^ I'm the complex conjugate amplitude. The direct terms contain singularities when 
partons ^ and j are collinear, and when parton j is soft but not necessarily collinear 
with parton £. Since our definition of the amplitude in the singular limits is based 
on the use of physical polarisations for the final-state partons (see Eq. f HTjl below), 
the interference terms do not exhibit purely collinear singularities. They do, however, 
contain soft and soft /collinear singularities. 

Let us first consider the direct terms , which are products of a splitting amplitude vi 
times a complex conjugate splitting amplitude v}. Summing over the daughter parton 
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spins and averaging over the mother parton spins, leads to the spin-averaged direct 
sphtting functions Wu which will form part of the subtraction terms 

= ^ ^ \M{pJ}m+i,Sj,se,Si)\'^ . (45) 

In d dimensions, when the mother parton is a gluon, the averaging factor becomes 
1/(2(1 — e)). An additional symmetry factor 1/2 has to be introduced for a final state 
g ^ g + g splitting. In the soft-gluon limit, Wu can easily be obtained from Eq. (jS]), 

TTTcikonal Pi- D{pj,Q) - pi 

Wu = 47ra, \ , 46 

[Pj-Per 

where flavour-dependent averaging factors are already taken into account. The trans- 
verse projection tensor D^j, in Eq. (H^ is given by 

n ^- A^ p';q'' + q^p'^ q'p';p', 

Pj ■ Q [Pj ■ Q) 

It will be convenient to define a dimensionless function F: 

F = ^Wu, (48) 

and we then have 

F-i, ^ " yy"'^"''^^ ^ P^ • DjPj^Q) -Pe ^ ■ Q _ Q^pe-Pj ^^g^ 
4:7!- as Pe-Pj Pj-Q {Pj-QY' 

For the direct terms, the colour part of the squared amplitude factorises 
(-M({p,/Wi)l-M({p,/}m+i)) 

~ C(/,, fMMi{{p, f}mm{p. /Wl)^/({P, f}m+l)\M,{{p, f}m)) (50) 

with overall colour factors 

= {Q,g),{9,<l) , 

CifiJj) = {Ca ifejj) = {g,g), (51) 

= (?,?)• 

We now turn to the interference terms, where parton j is emitted from parton i in 
the amplitude and from parton k ^ i in the complex conjugate amplitude (Fig. [2]). As 
mentioned above, the interference terms contain the soft and soft /coUinear singularities, 
and the splitting amplitudes are given by the eikonal approximation. The interference 
splitting function thus becomes 

Wek = t;,"({p, fU+i, se, se) °°^'({p, fU+u 4, s,)* 6s„s, 6s„s, • (52) 




There is an ambiguity in the momentum mapping to be used in Eq. fl52l) . One could 
associate Wek with the splitting of parton i and the momentum mapping for that 
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I 

\Mi{{p,f}m+l) > , <Mk{{p,f} m+1 J I 




I 
I 



Figure 2: Interference diagram: parton j is emitted from parton i in the scattering 
amplitude and parton j is emitted from parton k in the complex-conjugate scattering 
amplitude. 



splitting, or with the splitting of parton k and the corresponding momentum mapping, 
or one could use an average of both. To define the average, we introduce weight 
factors prj which redistribute the splitting function to the two possible mappings 



W, 



ik) 

Ik 



(53) 



Here W^,^ and W'^'^' denote splitting functions with the mapping corresponding to the 
splitting of parton £ and fc, respectively. Furthermore we have 



r(k) 



Atk{{p]m+l) + Akt{{p] 



(54) 



for any fixed momenta {p}m+i- The conceptually simplest choice Ai^ = A^i = 1/2 
has been adopted in Ref. [IQ]. However, we would like to set up our NLO scheme 
in a way that facilitates matching with a parton shower including quantum and spin 
interferences. Thus, as advocated Ref. |12], we will adopt weight functions which have 
more favourable properties when used in the formulation of parton showers. 

Summing over the two graphs with interference of gluons emitted from partons i 
and k we obtain a term 



W, 



IkM 



+ 



Aek 
Aki 



tk + Wkt t\ ® U 



4 + 

tk + wl:ftl®te 



(55) 



Let us consider the part proportional to A^k- The spin-averaged splitting function 
oc Aik can be simplified to 



1 

2 L 



tl te + tl tk 



W. 



ek 



where 



W 



ek 



4:7ias2Aik- 



i ■ D{pj,Q) -pk 
Pj-PePj-Pk 



(56) 
(57) 
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We now combine Wik with the direct sphtting function squared Wu, which has the 
same momentum mapping and which comes with a colour factor ® t^. Invariance of 
the matrix element under colour rotations implies 



4 ® = - ^ ^ [tl ® + tj ® tfc , (5 



and the complete contribution obeying one particular mapping is then given by 

1 



2 L 



tl (S) te + tl tk 



[Wee - Wek] ■ (59) 



Following [41j, we now add and subtract the soft-gluon approximation to the direct 
splitting function Wee according to 

Wee - Wek = {Wee - W^^^^^') + (w^f - Wek) , (60) 

^^^^~ClkOH3^1 ______ 

where VT^^ is given in Eq. f l46p . The first term in Eq. f lBU]) only has a collinear 
singularity, while the soft and the soft /collinear singularities are contained in the second 

term \W^n^°^^^ — Wek]- This second term can be written as 112] 



)2 



W^"^^' - Wek = 4™. A',, , (61) 



where Pek = Pj ■ PePk — Pj ■ PkPe, and A'^j^ is defined in |12] in terms of A^k- Several 
choices for A'^^ have been proposed in Ref . ; all results given here have been obtained 
using 

Pj -PkPi -Q + Pj -PiPk ■ Q 

specified in Eq. (7.12) of Ref. ^J. As argued in [12], the choice in Eq. f l^ has various 
favourable features when used in the formulation of a parton shower. 

The general form of the interference spin-averaged splitting function is then given 

by 

ATT7 r—e\kona\ jjj 2pe-pkPl-Q 

ATy = W^(^ - Wek = ^T^ois 7 — ^—r — rv • (63) 



Pi ■ Pj (pj -Pkpe-Q+Pe-Pjpk-Q 



The only singularity in Eq. fl63|) arises from the factor pe ■ pj in the denominator; the 
interference term is constructed such that it vanishes for pj ■ — )■ 0. Note that the 
interference term AW only needs to be considered if the emitted parton j is a gluon. 
If parton j is a quark or antiquark, this term vanishes. 



3.3 Subtraction terms 

The subtraction terms Ve in Eq. (|1|) are constructed from the splitting functions that 
describe the (m + l)-parton matrix-element squared in the soft and collinear limits, as 
derived in Section 13.21 Results for the splitting functions have already been presented 
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m Ref. Bill Integrating the subtraction terms over the one-parton unresolved phase 
space yields an infrared- and coUinear-singular contribution which needs to be combined 
with the virtual cross section to yield a finite NLO cross section. In this section we 
present results for the subtraction terms and their integrated counterparts Vi. We 
first consider the coUinear contributions proportional to W u — W^^^ "^'^ and then discuss 
the soft singularities contained in the interference terms proportional to W —Wek, 



see Eq. (|60|). 

In the formulae presented below, we leave out a common factor Anas in the ex- 
pressions for the squares vj of the splitting amplitudes; however, the expressions for 
the subtraction terms and the integrated subtraction terms Ve contain all factors, 
so that they can directly be used in the final formulae presented in Section 13.41 below. 
Some of the integrals have been calculated using the Mathematica package HypExp 
[191150]. 

3.3.1 Final-state collinear subtractions 

In this paper we consider processes with only up to two massless partons in the final 
state, so that ai = 1 {c.f. Eq. (|9])). We use the labelhng l^jjj. and Vjrjj. for a process 
with the splitting pe — )■ pe + Pj- For final-state collinear splittings, the subtraction 
terms can be expressed through the variables 

y = TT and z = , 64 

with 

Pe = Pe+ Pj, rii = {Q -Pi) = p^, Pe ■ Q = Pe ■ Q - Pe ■ Pj , (65) 

where Pm is the non-emitting final-state parton in the process. 

The integration of the subtraction terms over the one-parton unresolved phase space 
makes use of the phase-space factorisation 

d{pj}m+l 9{{pJ}m+l) = [d{pj}m] d^p9{{p, f}m+l) , (66) 

where g{{p, f}m+i) is an arbitrary function. The definition of the unresolved one- 
parton integration measure is 

d^, = dy e{y^^ < y < |/_) A'^'^^ ^ ^ 2 S^p'^) ^ 2 n 5+(^|) 

X (2 n)' 6^"^ (pe + Pj-\pe- o) ■ (67) 



^Note that there is an error in Eq. (2.38) of [U], which should read 



F = K 



cikoiial 
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Here y^in = for massless partons and Vmax is given by Eq. (fT2|) . For = 1, Eq. ([6 
reduces to 

= dyO iy(l - y)) (1 - 2.i*(p?) ^ 2. 

X (2 tt)'' (p, + ft - (1 - t,)p, - t/Q) . {ei 



Subtraction terms and integrals 
qqg, qqg 

The squared splitting amplitude for final state qqg couplings in the case of massless 
quarks is given by 



where 



'^qqg ^cik 



eik 



y 



4 



y{peQ) [l-y 

1 + y 



Feik + (l-e)z 



y 



y + z{l-y) {y + z{l-y)y 
Thus we have for the subtraction term 

11 Anas ^ / 2 „,2 \ 
The coUinear part of the integrated subtraction term is given by 



v; 



coll ^ 47ra, ^ 
q<ia 



27rfi^Y ( 1 



-. — 



4 7r r(i-£) V ViQ 



14 + - TT^ 



(69) 



(70) 



(71) 



(72) 



gqq, gqq 

The gqq splitting function for massless quarks is given by 

2 



am 



ypiQ 



{l-e-2z{l- z)). 



(73) 



Averaging over the helicity of the incoming particles gives an additional factor 1/(2 (1 — 
e)), and we obtain for the subtraction term 



911 



2{l-e) 

Integrating this over the unresolved one-parton phase-space, we obtain 
4:71 a. 



(74) 



V, 



,.2e. 



911 



2(1-8)' 





/27^/i2^ 


e 


1 8" 


TT T{l-e) 






~3e " 9 



(75) 
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The total (unaveraged) splitting amplitude squared is given by 
1 



999 



with 



and 



2 iPe ■ Pjf 

Pi ■ Dj ■ Pe 

Pj ■ De ■ Pj 
kl 

Tr [De ■ D,] 
A 



{{d-2) [pe ■ D, ■ Pi + pj ■ Di ■ pj] - klTi [De ■ D,] } , (76) 



y + z{l-y) 

'^ype-Q 

-2y z{l- z)pe- Q, 
d-2-2A + A^ 



y + z{l-y) - 



y 



y + z{l-y)_ 

y 

l-z{l-y) 



(Pi ■ Pj) 



2y 



(77) 



(7J 



(Pi ■ Q) iP, -Q) (?/ + ^ (1 - y)) (1 - 2; (1 - y)) ■ 

Instead of using this as a subtraction term, however, we proceed in a different way [H] : 
in order to well separate the singularities in the triple-gluon final state, we will use a 
slightly modified splitting function, where all soft divergences originating from particle 
i are transferred to the subtraction term where pi and pj are interchanged. This can 
be achieved by subtracting the term 



d-2 



where ^2,3 are defined corresponding to Eqs. (2.40)-(2.42) in [IT]. In the end, we obtain 



[pi ■ Dj ■ Pi - Pj ■ Di ■ Pj] 



(79) 



999 



7;2 _|_ n,2 
999 999,suh 



2 {pe ■ Pj] 



{2 {d-2) Pi- ■ Pi - klTi [Di ■ Dj] ] , (80) 



which is the unintegrated subtraction term for each gluon emission. The first part 
is the unaveraged eikonal splitting function; if we combine this with the interference 
term, we have 



~2 2 

999 eik 



'^{Pe-Pjf 

and the coUinear subtraction term reads 



^' Tr[A-D,] = [d-2-2A + A'] 



D 



coll 

999 



4 TT a. 



y {Pe Q) 



2(1-6) 

Integrating and taking all averaging factors into account, we obtain 

'27r/i2N^ 



(81) 



^2) 



v: 



coll 
999 



999 eik 



Qj r(l-e)2vr 



55 



^3) 
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Figure 3: Example for the initial state splitting q{Pa)(l{Pa)g{Pj), which has to be applied 
for processes with q + X — )■ Y in the m-particle and q + X — )■ g + Y in the (m + 1)- 
particle phase space. Our notation follows backward evolution, with the first flavour 
denoting the parton which participates in the hard interaction. 



3.3.2 Initial-state collinear subtractions 



In the following, the parent parton will always be pa, for processes with two initial-state 
partons, the corresponding formulae for the other incoming parton can be obtained from 
interchanging a 6. We label the processes according to backward evolution; e.g. gq q 
denotes a process where the gluon participates in the hard interaction. Therefore, 
this splitting function needs to be applied for g + X — )• F in the m-particle and 
q + X — )■ Y + q in the (m -|- l)-particle phase space. See Fig. [3] for an illustration for 
qqg splitting. 

We will use the following variables to describe in the initial-state subtraction terms: 



X 



1-2^ and y = 2P^-^^ 



4) 



g2 g2 

and we also define y' = y/(l — x). 

The momentum pj can be expressed in terms of these variables by introducing a 
Sudakov parametrisation 

pj = apa + f3pb -k±, (85) 

where 

1 — X — y 



a 



X 



and (3 = y. 



(86) 



From p| = it follows that k'j_ = —2 a Ppa - Pb- If we define the z— axis by the direction 
of the incoming beams, we can furthermore specify 

/ 

l-2v 
2^v{l- v) 



ki 



\ 







7) 



where v parametrises the additional angle in the interference terms; pj is then com- 
pletely parametrised in terms of x, y', v, 2pa ■ pt, and we can easily reconstruct 

1 , , 

-Pa +Pb- Pj (88) 



K 



X 
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required for the inverse transformation given in Eq. (!3T]) . 
The phase space factorises according to 

d{pj}m+l 9{{pJ}m+l) = [d{pj}m] d^p9{{pj}m+l) , 

with the 0?- dimensional integration measure 

A more exphcit form of the integration measure is given in Appendix [Bl 



Subtraction terms and integrals 
qqg, qqg 

The unaveraged qqg splitting function is given by 

vL„ - C;,, = z — (d - 2) . 



qqg eik 

Including all prefactors, we obtain 



xypa ■ Pb 



11 Anas ( 2 „,2 \ 



qqg 



and 



v; 



■coll 
999 



Cf j ^i^'d^p - ^'eik) 

4:71 fl"^ 



27r r(l-£) \2pa-pJ Jo X 



e 



- In X + 2 In (1 - 



gqq, gqq 

The gqq splitting function for massless quarks is given by 

,2 



(^-2) ^ A{l-{x + y)) 
V - = + 

yPa-Pb Pa-Pby{x + yy 



and we have 

Integrating, we obtain 

4 TT as 



V 



gqq 



V, 



999 



CF^^ I d^pV 



,2 

'J gqq 



1 a.. 



r(l-£) 27r \2pa-Pb 



Anas 2 
C^^%9 



2 \ e /•! 



11 / 1+ il-x) 
e X \ X 



+ 9{ 



with 



g[x,e 



0\ 



- 2 d - x) 



— In X + 



2 ln(l - x) 



x^ 



((1-xf + 1). 
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qgq, qgq 

The gqq splitting function for massless quarks is given by 



and we obtain 



+ 

xyPa-Pb xypa-Pb 



{{x + y)-l) 



cim 2 (1 - £ 
The integrated sphtting function reads 
4 TT a 



s rr 2 



V, 



191 



2(1 -el 



1 tts 



Til-e) 27r 



4 TTyU' 

'^Pa-Pb 



2 ^ ' dx 
^ 



(9^ 



(99) 



-- (x^ + (1 - xf) + 5'(x,£'^ 



(100) 



with 



g{x, £°) = (1 - x) (5 X - 1) + ( 2 ln(l - x) - In x) [x^ + (1 - xf] . (101) 



For initial states, the ggg splitting function is given by 



2 2 

999 eik 



2{l-x-y) 

y Pa-Pb 



(rf-2) 1 + 



1 



4(1 -x-y) 



x + yY) pa-pb{l-x){x + y) 



and 



coll 



47r a. 



999 2 (1 - e 
After the integration, we obtain 

4 vr a 



■Ca ( 



999 ^eik 



J- 



(102) 
(103) 



v: 



coll 
999 



^— ^ /i^^ J dip {vlgg - t;cik) 



a. 



X 



rfx 2 (1 - x) 



1 + x^ 



x^ 



+ 2 ln(l -X 



r(l - e) 27r 

X (2 — x) 



(l-x)2 



4 vr /i 
2 Pa -Pft 

In X 



2 \ £ 



+ 



X 



1 — X X^ 



(104) 



3.3.3 Soft and soft/collinear subtractions 

We now turn to the discussion of the soft and soft/collinear singularities which arise 
when the emitted parton is a gluon and which are contained in IF^^ — Wik, as given 
in Eq. fl63|l . There are contributions from initial-initial-, initial-final-, final-initial-, 
and final-final-state interference terms, which we now discuss in turn. 
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Initial—initial-state interference 



In terms of the variables x and y defined in Section 13. 3. 2^ the initial-initial-state inter- 
ference term is given by 



1 



4:71 a. 



-AW, 



ab 



I -x)ypa ■ Pb 



(105) 



so that 

Vi = QAWab. (106) 

For the simple scattering processes we consider in this paper, the colour algebra fac- 
torises and leads to the colour factors Ci = Cp {Ca) for the qqg (ggg) splitting func- 
tions. 

Integrating this over the phase space of the unresolved parton, we obtain 



Atc fi"^ \^ as 1 



2pa-pbJ 277 r(l - e] 



1 

X I dx 

'0 



-7 ---6{l-x)-2j^ 



+ 4 



ln(l — x) 
1 — X 



(107) 



Initial final-state interference 

For the initial-final-state interference, we obtain 



-AW, 



ak 



2 {Pa ■ Pk) 



{Pa ■ Pj) {Pj ■ Pk)s + 2 {pa- Pj) {pk ■ Q) 



and 



V^l, = C AW, 



ak • 



Using the initial-state integration measure given by Eq. fl90|) . we have 



1 



4:71 fi"^ \^ as 

2pa ■ Pb ) 2^r(l-£) 



1 



X / dx < — 6(1 — x) , 



+ 5(1 — x) In^o 



+ 4 



ln(l — x) 
1 — X 



2 In X 2 
— + 4(5(1 - x) (In So) (In 2) + - hn{x, z) 

, i — 3^)4- TT 



(108) 



(109) 



(110) 
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with 



Jfin(a;, z) = Tc (5(1 — x) X 



j;2 



In 



2y' +~Zo-2y'~Zo + ^Ay'^ (1 - ~Zo) + 5, 



i:2 



In 4 In 5o + ^ (Li2 (1 - ^o) + In^ ^o) 



+ 



1 



(l-a;)+ Jq y' 



dv 



^Jv{l - v) N{x,y',z,v) 



- 1 



:iiii 



and 



Pa-Pk 

xpk-Q' 



Pa -Pk 
Pk-Q 



N 



Pj ■ Pk 



Pk-Q '^-x 



+ y' 



(112) 



Note that z imphcitly depends on x, y and v through momentum mapping. Further- 
more, Pk has to be obtained from pk using the transformation specified in Eq. (13T|) and 
with K given by Eq. (188|) : pj is given by Eqs. (185|) . (186|) . and (IHTl) . 5 = corresponds 
to a singularity in the m-particle phase space; this singularity should be excluded by 
an appropriate infrared-safe jet function. As before, Ci = Cp (Ca) in the qqg (ggg) 
splitting function. 



Final— initial-state and final— final-state interference 

The final-initial and final-final-state interference terms have the same structure, with 
the only difference of pk = pe (final-state particle) and pk = Pa (initial-state particle) 
for the initial- and final-state integrals, respectively. 

For interference terms with a final-state emitter and final-state momentum map- 
ping, we obtain 



1 



47ra 



-Aiy 



2 {Pt ■ Pk) {pe ■ Q) 



{pi ■ Pj) ( {Pj ■ Pk) {pi ■ Q) + {pi ■ Pj){pk ■ Q) 



and 



V'^ = a AW . 



(113) 



(114) 



These are the general expressions for interference terms using the final-state mappings; 
if the spectator is an initial-state particle, pk = Pa- 

We obtain for the integrated interference term with final-state spectators 



n 2 \ £ 
2 jj, TT \ «, 



^Pe-Qj TT T{l-e 
with Ci = Cp {Ca) in the qqg (ggg) splitting function. 



1 1 TT^ 

+ T + 3 , 



2e 



(115) 
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For interference terms where the spectator is the initial-state parton pa, we obtain 
for the integrated subtraction term 



2 /i 77 \ a 



Pe-Q J IT T{1 -e] 



a 



X 



1 1 
+ - 



2e 



TC 



1 + - In (fio + 1) 



+ 3 - 2 In 2 ln(ao + 1) + 

6 TT L 



4?(«o) + , (116) 



with 



1 



^ du , , 
— <^21n2+ ^ 
u I x/l + 4ao(l + ao)M2 



X In 



[l-u) 



{l + 2aoU + a/1 + 4ao (1 + ao) 



+ 2 In 2 ln(l + ao) 



hr (1 + ao) 



5.. / ao 

-Li2 

2 V «o + 1 



Lio 

2 



ao 



cto + 1 



TT 



X 



x{l — X + ux [{1 — ux) a + 2]) 
k{u, X, a) 



l + AaoU^{l + ao) 



We have introduced 



k'^{x, u, a) = [(1 + ux — x){z — z) + ux ((1 — ux) a + 1)]^ 
+ 4:ux z {1 — z) {1 + ux — x) ((1 — Ma;)a + 1) 



and 



X {1 — u) I 



Pa ■ ne 



1 — ux 



Pa-Pt + yPa ■ ne 



ao = d{y = 0) 



Pa ■ ne 
Pa- Pi 



(117) 



:il8) 



(119) 



(120) 



Note that the treatment of interference terms significantly differs from [31j; here, our 
choice of momentum mapping leads to more complicated integrated interference terms, 
the finite parts of which we choose to evaluate numerically. 



3.4 Final expressions 

Let us finally collect the formulae that are needed to address the scattering processes 
considered in this paper. The NLO parton level cross section for any collider process 
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is given by the sum of a^^^ and cr^f,^", with 



,NLO 



(J. 



ah 



dcrah{Pa,Ph) 



a. 



NLO 
ah 



m+1 



d(TabiPa,Pb) + / da^i,{pa,Pb) + / da^^iPa, Pb, ^l) ■ (121) 

J m Jm 



The hadronic cross section is obtained from the parton level cross section by convoluting 
with parton distribution functions. The collinear counterterms J,^da^i,{pa,Ph, fJ'^?) are 
needed to absorb initial-state collinear singularities into a re-definition of the parton- 
distribution functions; in the MS scheme, they are given by 



/ da'^^{pa,Pb,^il) = ^TTTT — T ! dx f da%{xpa,pb) 
Jm 27r r(l -e) Jo Jm 



1 / 47r/i 



2\ e 



+ 



1 /47r^2x- 



(122) 



Here the P°'^{x) are the Altarelli-Parisi kernels in four dimensions |15]; their explicit 
form is given in Appendix O We then obtain for the parton-level NLO contribution 



(^ab'^iPa, Pb, = / [daj^iPa, Pb) " da^i^{pa, pb) 

Jm+1 

d(yIbiPa,Pb) + J^d(J^f^{pa,ph) + da^biPa,Pb,fJ'l) 

where da^f^ + da^^, can be written as 



+ 



(123) 



e=0 



dc^abiPa, Pb) + da^biPa, Pb, fip) 



d(Tab{Pa,Pb) ® I{e) + dx da^b{xpa,Pb) ® [K^ixpa] + P{x,fil) 

Jo Jm 

+ dx daj^^{pa,xpb)^[K\xpb) + P{ ) • (124) 

^0 Jm 

This equation defines the insertion operators I{e), K{x), P{x; fxp) at the cross sec- 
tion level, where we follow the standard notation introduced in Ref. [31]. Eq. fll24p 
can be divided into two parts: the first part is the universal insertion operator I{e), 
which contains the complete singularity structure of the virtual contribution and has 
LO kinematics. The second part consists of the finite pieces that are left over after 
absorbing the initial-state collinear singularities into a redefinition of the parton dis- 
tribution functions at NLO. It involves an additional one dimensional integration over 
the momentum fraction x of an incoming parton with the LO cross sections. 

Obviously, the cross sections and all observables have to be defined in an infrared- 
safe way through the introduction of jet functions; one should thus replace a^^ and 
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i^NLO QYoss sections 

Cr^^ = J dPSm{Pl, ■■■ ,Pm.) \Mm{Pl, " " " ,Pm)f Ff'\pi, ■■■ , Pm) 
^NLO ^ J dPSm+liPl, ■ ■ ■ ,Pm+l) \Mm+l{Pl, " " " ,Pm+l)f F^"'^^\pi, ■ ■ ■ , pm+l) 

+ j dPSmiPl, ■■■ ,Pm) \MmiPl, ' ' ' , Pm) iLc-loop Pj^\Pl^ ' ' ' , Pm) ■ (125) 

In general, the jet function may contain ^-functions (which define cuts and correspond- 
ing cross sections) and 5- functions (which define differential cross sections). Infrared 
safety now requires that 

pbn+l) [p^^... ^p. = Xq,--- , Pm+l) i^j""^ (Pl, ■ ■ ■ , Pm+l) H A ^ 

F!r^^\pi, ..,Pi, ..,Pj, ..,Pm+l) -> F^"'\pi, ■■,P, ■■,Pm+l) if Pi Zp, Pj ^ (1 - Z)p 

Ff"^ (pi, ■ ■ ■ , p™) -> if Pi- pj ^ (126) 

The first two conditions of Eq. (11261) define the essential property of the jet function 
that the jet observable has to be infrared and coUinear safe for any number m of partons 
in the final state, i.e. to any order in QCD perturbation theory. The last condition of 
Eq. fll26p guarantees that the Born-level cross section is well defined. To summarise, 
we require that 

p{m+l) ^ p{m) (^27) 

in the singular limits. 

For processes where one or both incoming particles are leptons, the collinear coun- 
terterms are set to zero and the parton distribution functions are replaced by 6- 
distributions, i.e. /^y" = 5(1 — rji)^ 

In the following, we discuss the specific form of da^i^{pa,Pb) which corresponds 
to the subtraction term in the real-emission contribution of the process, as well as 
the integrated d-dimensional counterterm j^da^^^{pa,Pb)- In general, the subtraction 
term can be split in maximally four contributions for processes with maximally two 
final-state particles in the leading-order contribution; in our scheme, each of these 
contributions requires exactly one momentum mapping. We then have 

dcTabiPa, Pb) = j,'"(Pa, Pb) + da^biPa, Pb) + da^^^ {pa, pb) + da^^^' {pa, pb) , (128) 

where ki^2 now label the momenta of the outgoing particles in m-particle phase space. 
The associated real-emission contributions can be obtained from considering possible 
splittings in flavour space: the sum over flavours in each of the above contributions is 
such that (/a, /fcj = {fa + fkj,fki) all combinations, i.e. we only consider flavour 
mapping where fa = fa + fkj is physically allowed. We will denote this by a delta 
function in flavour-space mapping S^.^j. From initial-state splittings, we then have for 
the single contribution with fixed initial and final-state flavours (in the following, we 



^This of course only holds for the discussion of higher-order effects from strong interactions; for 
electroweak processes, the inclusion of structure functions can not be neglected; c.f. 51j and references 
therein. 
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omit the jet functions for notational reasons; however, full expressions should always 
be read according to Eq. fll25p where all jet functions are included): 



d(r^%Pa,Pb) = \ l'^9'lM)^g■,g,g^+'^999(.P^)^9■,9,9^]\^Born,g\\xpa,Pb■,Pn) 

+ [T>ggg{Pi)Sg.g^g^ + 'Dggg{Pi)6g.g^gXPi)] \M BoTuj^ P a , Pb] Pu)^ , (l29) 

where Nm+i incorporates all symmetry factors of the m + 1 process and $m+i = 2 s is 
the respective flux factor. D{pi) now signifies that pj = pi, and A^Born,g(g) corresponds 
to the underlying Born matrix element for the process Pa + Pb Yin P^ with an 
incoming gluon (quark) such that fa = 9 (q) and all other flavours unchanged, and 
where n ^ i labels the momenta of the final state particles. The momentum mapping 
is here given by Eqs. ( 136|) and ( 138|) respectively. Note that, while Vggg and Vggg 
contain coUinear singularities only, Vggg and Vggg are split into a coUinear and soft and 
interference term as discussed in Section 13.3.31 

-I^Up^) = + ^'^P^^Pb)^L9 + (130) 

where D{pi,pk) now denotes an interference contribution where pk acts as a spectator. 
An equivalent expression holds for Vggg{pi). The second term sums over all final-state 
particles with k ^ i which are gluons. da^^{pa,Pb) is obtained from da^^"'{pa,Pb) by 
interchanging a ^ 6 in the corresponding expressions as well as in the flavour mapping 
functions. All subtraction terms given here refer to the expressions in Section I3.3.2[ 
and to the expressions for initial-initial interference terms for V^\pi,pb) and initial- final 
interference terms for V^\pi,pk) in Section [3.3.31 

For the final-state splittings, we have a similar expression 



d(T^''\Pa,Pb) = V { [Dggg{pi)Sg.g,g^+Vggg{pi)5g.g^g^]\MB0Tnj'^{Pa,Pb;Pn) 



+ PgqgiPi)^q;g,Q^ + '^QqgiPi)^Q;q,gAPi)] |-^Born,gP(Pa, Pn) | , (131) 

where now the delta functions in flavour space are defined as Sk^-p^,pj, and A^Born,g(g) 
corresponds to the underlying Born matrix element for the process Pa + Pb Pn+Pki, 
where now f^. = g (q) and all other flavours unchanged, and where n ^ ki labels the 
momentum of the remaining final-state particle. As before, T>qgg and T>ggg contain both 
coUinear and interference terms, and we have 

-DgggiP.) = ^^^(p.) + T^^" iP^^ P^^)^ Lg + ^ ^^^^^^ ) ^/-« (1^2) 

k=a,b k^i 

and the same for Vggg. The mapping needed to define the subtraction terms is given by 
Eqs. f l2ip and f l24p . All subtraction terms given here refer to the expressions in Section 
I3.3.1[ and the expressions for the final-initial interference terms for V^\pi,pa^b) or the 
final-final interference terms for V^\pi,pk) in Section [3.3.31 
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The integrated counterterms are generically given as 

/ da^h{pa,pb) = / da^,;^{pa,pb)+ / da^\pa,pb)+ / da^'^^pa, Pb) + / da^^\pa,Pb) ■ 
Jl Jl Jl Jl Jl 

(133) 

The collection of the integrated counterterms is then straightforward: for each term 
which has been subtracted in the real-emission part, the corresponding integrated con- 
tribution to J, K, P needs to be added to the virtual contribution as in Eq. f ll23p . 
Finally, note that the expressions given in the following sections are derived on the 
matrix-element level: 

^ ^ jv\M\l = V\Mt; (134) 

on the cross-section level, we additionally have to take the flux and combinatoric factors 
into account, such that 

j^dat^iAPa.Pb) = ^ jviMt = ^V\M\l, 

[ [ dat^,,, = iV^+x / ^V\M\l = %^xV / da^. (135) 

Jm Jl Jm -''m J m 

Here the factors A^m, Nm+i account for possible symmetry factors of the specific process. 
We then have the relation 

Y.^ = \P^{I + K + P) (136) 

m+l 

between the integrated splitting functions V derived in the previous section and the 
insertion operators /, P as defined in Eq. 01241) . 



4 Application to physical processes 

In the last section, we have presented the subtraction terms for the real emission 
and their integrated counterparts for processes with up to two particles in the final 
state. We will now use these expressions in well-known processes at NLO, showing that 
indeed the singularity structures of the real-emission terms and one-loop contributions 
is reproduced by the squared and averaged (one-particle integrated) splitting functions 
of the parton shower [ID] in the singular limits. We validate our scheme by showing 
that the application of our subtraction terms reproduces standard results from the 
literature for all processes. In more detail, we discuss single-W production (initial-state 
qqg and qgq splittings, initial- initial state interference term), dijet production at lepton 
colliders {qqg final-state splittings, final-final state interference term), gluon- induced 
Higgs production in an effective theory description {ggg and gqq initial-state splitting 
functions, initial-initial state interference terms), Higgs decay to two gluons {ggg and 
gqq final-state splitting functions, final-final state interference terms) and a subprocess 
of deep-inelastic scattering (DIS) {qqg initial and final-state splitting functions, initial- 
final and final- initial state interference terms), c.f. Table [H We find that the splitting 
functions for the parton shower as described in jlHl SH 112] , in combination with the 
momentum mapping, are well suited as subtraction terms in NLO-QCD calculations. 
The numerical evaluations of phase-space integrals in this section have been obtained 
using routines from the Cuba library |52] . 
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Splitting function 


Process 


gqq, initial state 


Higgs production 




sinde W DIS 


qgq, initial state 


single W 


llllljlCXl OuCXuv!' 


l-Tio"o"c; nrnnnpfinn 




c;Tno"lp vV rT 1 (Tcrc "nTonnp't'io'n 


liilllcLl iilicLl liiLci ici ciiL-c 




gqq, final state 


Higgs decay 


qqg, final state 


Dijet, DIS 


ggg, final state 


Higgs decay 


final final interference 


Higgs decay, Dijet 


final initial interference 


DIS 



Table 1: List of all splitting functions presented in Section [3731 and test processes used 
for the scheme validation in Section |H 

4.1 Single- W production 

We start with a simple process: single-W production at a hadron collider. The tree- 
level process here is given by 

(1{Pi)q{P2) — > w, 
and real-emission processes include both quark- and gluon-induced cases 

(1{Pi)q{P2) — > Wg{p3), g{pi)q{p2) — > W q{p3). 

This process contains both the qqg as well as gqq initial state splittings given by 
Eqs. (192|) . (199|) as well as the initial state interference term Eq. (11061) . All results for 
tree-level and real-emission matrix elements are well known and have been taken from 
the literature (e.g. [S3])- 

Tree-level contribution 

The squared matrix element for the tree-level process qq' W is 

I Mb r= I V,,, r Ml. 

where we have averaged over initial-state particle spins and colours. The one-particle 
phase space is given by 

J dPSi = 27r5+(s-M^) 
with s being the partonic center-of-mass energy. 
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Virtual correction 

The virtual contribution in the MS renormalisation scheme is 



2Re (MbM:,,,) = \Mv^' 



' virt / 

Mb \' ^ 



27r^"r(l -e) \ Q 



2\ e 



2 3 



e 



8 + 7r2 + C(£) \ ,(137) 



and we have 



„virt 



1 / ATljJ, 



27r r(l -e) V Q 



,2\ ^ 



2 3 



+ vr > cr . 



(138) 



Real emission 

The matrix element for the quark-induced NLO real-emission process qq' — )> Wg is 
given by 



Ml? 



V2GfM^ I V,,> P — 



32 £2 + -^2 + 2M^s 



9 



9 tu 



tu 



(139) 



where (7 and are related to the Fermi coupling constant Gi? by ^ = g^^, and 

we use s = {pi + ^2)^, i = iPi — PzY s + i + u = m^. The matrix element for the 
gluon-induced process can be obtained by crossing symmetry 



Ml? 



The two-particle phase space is 



-g T^as I Vqq' 



+ + 2M^i 



-su 



dPSi 



(27r)32p0^ (27r)32p!] 



PW -P3) 



Subtraction terms 

In the quark-induced case, we need to consider two subtraction terms corresponding 
to the gluon emission from q and q, respectively, which we label and ■ Their 
definition is given in Eqs. ( 192|) and Eq.f ll06p . We have 



4 vr a. 



Au 



+ 



i us 



IQsu^ 



tM'^ ■ t{u + tf {t^ + u^){t + uy 



(140) 



where we consider the case where pi = Pi, J-e. gluon emission from the incoming quark. 
The subtraction term can be obtained from Eq. (11401) by the replacement t ^ u. 
As the Born matrix element is constant and we use a unit jet function, no momentum 
mapping is required. 
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After subtraction of the counterterms, the final expression for the two-particle cross 
section at parton level is given by 



NLO{2} _ 



^ J dPS2{\M'^\' -{VI + VI)\Mb\'} = 0. (141) 
The splitting function in the gluon-induced case Vg is given by Eq. (199|) . and we obtain 

= Ysj dPS,{\M'^\'-Vg\MB?] 
Integrated subtraction terms 

In the m-particle phase-space contribution, we now have to consider the contributions 
from both quark- and gluon-induced real emission processes; note, however, that the 
poles from the virtual contribution are completely cancelled by the integrated qqg 
splitting function. In this case, the corresponding integrated subtraction terms are 
obtained from Eqs. (nUTj) . 

The collinear singularity needs to be cancelled by the universal collinear counter 
term as explained in Section 13.41 We then obtain for the quark-induced case 

/^C°^^^= /^^fj2V,,,) + ^rfaf, 

= / dal^i,{pa,Pb) ® lq{e) + dx da^f,{xpa,Pb) ® [Kg{xpa) + Pq{x,IJ,p)] 
Ji Jo Jl 

(143) 

The corresponding I, K and P terms are 



-(1 - x) Inx + 2(1 - x) ln(l - x) + 4x 



2tt r(i -£) 

2a; In a; /l + x^\ / 47ru^ 

In 



ln(l — x) 
1 — X 



il-x)^ \l-xj^ \2xpa-Pb 

We immediately see that the singularities in | A^y P and 1(e) cancel. 

The integrated subtraction term for the gluon-induced case is given by Eq. fllOOl) . 
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and combined with the corresponding coUinear counterterm, we have 



1 

Jl 



(144) 

dx I da^h{xpa,Pb) ^ [Kg{xpa) + Pg{x,fip)] 



where 



Kgixpa) = g^ ^,/ , I TniGx - 5x^ - 1) + [2 ln(l - x) - ln(x)] P^^ix) 
Zir L [I — e) I- 



2xpa ■ Pb 



P^\x) = Tn[x' + {l-xy], Tr 



1 

2 



Results 



The complete partonic cross section for the single-!^ production at NLO can then be 
obtained by combining Eqs. ffT38|) . ffTiHj) . ffMij) . ffTiT]) . and f lTi2|) : 

^NLO ^ ^virt ^ ^NLO{l} ^ ^NLO{l} ^ ^NLO{2} ^ ^NLO{2} _ 

Comparison with a calculation in the Catani-Seymour scheme reveals that in the two 
schemes the finite terms are attributed to diff'erent parts, but the complete contribu- 
tions of course agree. For a comparison, we implemented both schemes in a private 
code; Fig. O shows the relative difference between the two schemes as a function of the 
hadronic center of mass energy, convoluted with parton distribution functions. The 
schemes are equivalent, with agreement which is consistent with zero on sub-per-mil 
level. 
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0.35 



0.25 




0.15 



Figure 4: Relative NLO correction (o"nlo ~" o"Born) / o'Eom to single-W production at the 
LHC as a function of the hadron center of mass energy. The resuh was obtained using 
the CTEQ6M parton distribution function |54j . 



0.001 



5x10- 



-5x10-"- 



-0.001 




;[TeV] 



Figure 5: Relative difference between NLO corrections to single-W production using 
the Catani-Seymour (CS) and our subtraction scheme based on Nagy-Soper splitting 
functions (NS), as a function of the hadronic center of mass energy. The results agree 
at the sub-per-mil level, shown are the numerical integration errors. 
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4.2 Dijet production 



We next consider dijet production at NLO. In this process, the subtraction terms are 
final state qqg sphtting and the final-final interference terms, given by Eqs. fl7ip . fill 41) . 
The integrated counterparts can be obtained from Eqs. fl72l) . (11151) . We here closely 
follow the procedure in [31], i.e. we calculate the prefactor which relates the NLO and 
LO cross section via 

4 TT 

We average over the event orientation in the leading-order process; in this case, the 
momentum dependence of the Born contribution vanishes. 

Tree-level contribution 

We normalise the leading-order matrix element such that the phase space is given by 

2pi ■ P2 



J dPS, = j dy6{l-y), y 

a^o = \M2? I dy6{l-y)Ff 



J 

and pi,P2 are the momenta of the final state partons. 
Virtual contribution 

The one-loop matrix element in the MS renormalisation scheme is given by 
2Re(A^2-M:i,J ^ \Mv\' 

= \M2\^^Cf^, {-^ 8 + n^ + 0{6) 



2ti 'V{l-e)\ ) \ e 



= S^^-FTT^ - 7 - 8 + da-- . (146) 



and we have 

'''' " 27r"^r(l -e) \ ) \ e 
Real emission 

The matrix element for the NLO real-emission process 

e+e" q{pi) q{p2) giPs) 

in four dimensions is given by 



Q2 (1_3,^)(1_3,2) 

where we have used Xj = and | M.2 P is the LO matrix element for the process 
e~^e~ q + q. For the phase-space integration, we use 

/ dPS'i = — / dxidx2Q(xi + X2 - I). 
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Subtraction terms 



The calculation of the subtracted splitting function contains two contributions, Vi and 
V2, which correspond to emission from the quark and the anti-quark, respectively. 
They are given by Eqs. fl69|) and f lll4p . We obtain 



4 

+ 2 



Xi 

X3 



1 — X2\ 1 — Xi 



1 - X2 



1 

1 - X2 J (1 - Xi)Xi + (1 - X2)X2 



Xi 



(147) 



where we have used Q = pi +P2 +P3 = Q, X1+X2 + X3 = 2, and = is the square 
of the center-of-mass energy. The contribution V2 can be obtained from Eq. fll47p by 
the replacement Xi ^ X2- The final expression for the three-parton cross section is 
then given by 



a 



NLO{3} 



(^Cp)\M2 P f' dXidX2Q{Xi + X2-l)l—^^r^ tF 



(3) 
J 



+ 2 



I- X2 

Xi 



I-X3 



1 - X2/ (1 - Xi)Xi + (1 - X2)X2 



Ff (1) 



Xi 



X3 

1 — Xi J (1 — Xi)xi + (1 — a;2)a;2 



1 — Xi 

X2 



Ff\2) 



(148) 



where F) [i) signifies that pi = pi in the corresponding mappings in the jet functions. 
For any infrared-safe observable, Fj^\i) — )■ Fj^^ as Xj approaches 1, and the above 
expression is finite. For unit jet functions, it reduces to 



^NLO{3} ^ ^cJ---n'\a 
27r V 2 3 



LO 



(149) 
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Integrated subtraction terms 



The integrated collinear and soft subtractions are given by Eqs. fl72p and flllSp : 



as^ 1 



2\ £ 



Air a. 



■^oik ^soft^ 



47r/i 

27r^"r(l -£) \ Q 



1 



2\ £ 



(150) 



Combining these two contributions with the virtual cross section, we obtain a finite 
expression for the two-parton cross section: 



cr 



NLO {2} 



da^ + / da' 



e=0 



j dPS2{\ Mv ? +2 (V,,, - I M2 ?] F 

M2 ry dy6il-y)F 



(2) 
J 



(2) 
J 



(151) 



Result 

Summing Eq. fll49p and Eq. (115 ip . we obtain for unit jet functions 



cr 



NLO 



^NLO{2}^^NLO{3}^ 3«,^^^ 

4 TT 



LO 



which agrees with the standard result in the literature. Note that, although the 
angularly-averaged tree-level matrix element is independent of the parton momenta, 
non-unit jet functions need to account for the mappings from m -|- 1 to m phase space 
as specified in Section I3.1.1I 



4.3 Gluon- induced Higgs production 

In this section, we discuss the NLO corrections to gluon-induced Higgs production, 
where the leading-order contribution has been derived in an effective theory approach 
[55| [56| \57\ [58] . This involves ggg and gqq initial-state subtraction terms given by 
Eqs. (llOSp . (p5|l as well as the initial- initial interference term Eq. (I106P and their 
integrated counterparts Eqs. f lTUI]) . (1^ . ffTOTl) . 

The real-emission contributions in this process are 

99 9H, qg qH, qq gH, 
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such that the full NLO QCD cross section includes 



J d(Tqq^gH + Jdaqg^qH + J dagg^gH + J day + J dacon + J dach + Jdaes, 

(152) 

where J dach, J dacs are contributions from charge renormalisation and the matching 
of the effective to the full theory, respectively. In the following, we use unit jet functions 
throughout the whole calculation. 

Tree-level contribution 



The lowest-order cross section was first presented in ED]- From the QCD La- 
grangian, we obtain for the one-loop induced Hgg coupling the following expression 
for the total cross section, where the 5-function accounts for the one-particle final state: 



\A\ 



where 



rq)firq)) 



4M2 



Gp 

71 



(246 GeV)^ 



f{rq 



sm 



In 



l-r„ 



— it: 



if rq>l 



if Tq<l 



and s is the partonic center-of-mass energy, Mq the pole mass of the heavy quark (which 
we assume to be the top quark), and Mh the Higgs boson mass. In the limit that the 
top-quark mass is infinitely large, — )■ oo, A — | and 



(153) 
4 - 2e 



Note that the cross section is here given in = 4 dimensions; the d 
dimensional cross section is related by ct^q^ = ao/{l — e). 

The cross section to 0{al) for gg ^ H m the limit Mt ^ Mh can be obtained 
from the effective Lagrangian 



C 



eff 



a. 



where is the gluon field strength tensor and the index A is the colour degree of 



freedom of the gluon field (which runs over 1 ■ ■ ■ 8). All processes considered here and 
the corresponding Feynman rules have been obtained using this effective Lagrangian. 
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Virtual matrix element, charge renormalisation, and effective Lagrangian 
correction 

The one-loop matrix element in c? = 4 — 2e dimensions in the MS renormalisation 
scheme is given by 



2Re(7WLo-M:irt) = \M 



|2 



V 



\M^of ^^^(^)' r(l + e) (-1 + + 0{e) ) x 2 , (154) 



2 o^lMj, 



where 

5767r2t;2(i_e) 

is the (i- dimensional leading-order squared matrix element. For the total cross section 



avirt = ^0 — \ 1 r (1 + e) 



11 2 . 

-— 1 + -TT^ 



where now ctq denotes the four- dimensional LO cross section as given in Eq. f ll53p . We 
also need to take charge renormalisation into account; the charge counterterm in the 
MS renormalisation scheme is (see e.g. |M] 

= {^Zg) alf^ (156) 

where 

2 \ £ 



and Uf is the number of light quarks. Finally, we have to add the 0{as) correction to 
the matching coefficient which relates the effective Lagrangian to the full theory, 

11 

C^eff = ctq ■ (158) 

TT Z 

Real emission 

The real-emission contributions relevant for this process are given by 

9{Pa)9{Pb) ^ 9{pi)H, q{Pa)9{Pb) q{pi)H, q{Pa)q{Pb) -> 9{Pi)H, 
and the Mandelstam variables as defined as 

S = {Pa+Pbf , i = {Pa-Plf , U = Mjj-S-i. 

Pole cancellations occur between the purely gluon-induced real emission process and 
the virtual contribution; all other processes contain at most coUinear divergences which 
are eliminated by the addition of collinear counterterms. 
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For the real-emission matrix elements, we obtain for NLO gg gH in four dimen- 
sions 

\M{gg^gH)\ 



,2 al 32 Ml + s^ + + 



37r stu 

which has singularities when t — )■ or -u 0. The spin and colour averages yield an 
additional factor 1/2/(1 -e) x 1/2/(1 - x 1/8 x 1/8 = 1/256/(1 - e)^. 
The quark-gluon induced contribution is given by 

I A^/ rrM2 16 al S^ + V? 

\M{qg^qH)\ =- — 

which has a singularity when t — )■ 0. The spin and colour averages yield an additional 
factor 1/2 x 1/2/(1 - e) x 1/3 x 1/8 = 1/96/(1 - e). 

The purely quark-induced matrix element can be obtained from this by crossing 
symmetry, and we have 

\M{qq^gH)\ 2 — , 

9 irv^ s 

which is completely finite. The spin and colour averages yield an additional factor 
1/2 X 1/2 X 1/3 X 1/3 = 1/36. 

Subtraction terms 

We only need to consider subtractions for the gg and qg induced cases; these are given 
by Eqs. f ll03p . fl95|) . and f ll06p respectively. Combining them with the real emission 
matrix elements and integrating over phase space, we obtain 

a'^^o^^gg ^ gH) = J dPS^ {\M{gg ^ gH)\' - {D,,^ + D'^,) \Mlo\'} 



384s \ ti^v^ 



\ {Ml - sf ""KMIJ ^ 3i ^ Ml-s 



a^'^'^'^qg ^ qH) = 1- J dPS2 {\Miqg ^ qH)\' - V,,, \Mlo\'} 



(159) 



1728s VvrV J s J V^" ^ 

(160) 

As the Born contribution is constant and we are using unit jet functions, we do not 
need to apply any momentum mappings. 

The finite contribution from the quark-induced case is 



(161) 
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Integrated subtraction terms 

The integrated subtraction terms for the gluon- and quark-induced cases are given in 
Eqs. f ll04p . f l96p . fll07p . After addition of the collinear counterterm, we have for the 
gluon-induced case 

Jo 

+ / dxdao{pa, XPb) [VgggiPb) + VabiPb)] + da'^ {pa, Pb, Hp) + dach 

Jo 

= dao{pa,Pb) ® lgg{e) + / dxdao{xpa,Pb) ® [K'^ggixpa) + Pgg{x,IJ.p)] 

Jo 

+ (a ^ b) (162) 
where V{pi) signifies that pi = pi in the respective integrated counterterm and where 

a, 1 ( A'Kii'^y (2 Ti' 



2ti T{l-e) \ Q2 
a 1 



-2 — + x(l-x) + In' ^ 



2^f 



X J \2xpa-pb 

_L_5(i_^) (^}jcA-lnfTn^ ln(4vr), 



Pgg{x,fll) = ^ T^,/ Pgg{x)\n^^'^^ 



with 



Agg{x,e^) = ^^(-J—^^ -2x(l-x)lnx + 4(l-x)ln(l-x) (^i^^^ 



+2 X 



+ 

1 — X 

X 



Pgg{x) = 2C^ + + (^yC^-^n/Tfi 
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Integration over phase space yields 

.3 \ I 



a 



NLO,l 
99 



576 [v^Ti^) '^^r(l-£) 



X 



4 TTyU 



2\£ 



11 TT^ 

+ - + 1 -— 1 6{l-z)+4: 



e 



6 



ln(l - z) 
1-z 



+ 



-z {2-z{l-z)) In(l-z) 



+2 [z^ -{1-z)] - 2 



{1-zy 



+ z''{l-z) + {l-z) 



In^ 



+ 2 



(1-^) 



+ (1-^) 



lnz+^J{l-z) (^^C^-^UfTn^ In^j, 



(163) 



where we introduced z — Mfj/s. For the contribution from the ggi- induced case, we 
obtain 



UOgg 



dxda^^{xpa,Ph) \ygqq{Pa) + '^'IhiPa)] + (^abiPa, Pb, 1^1) 

dxda^^{xpa,Pb) ® [Kgg{xpa) + P<j3(a;,Ai^)] , 



a. 



where 

^qg{xPa) 

P{x,i^l) 

P^X) = Cp 

Aqg{x,e^) is given by 



L_^{^„(.,e»)C.-n,Wl.(^) 



27rr(l 

27rr(l-£)^^^^^^^''lv //I 
l + (l-x)2 



Aqg{x,e^) = ^ ^ - j;lnx + 21n(l -x) ' ^ ' 



X 



X 



Integration over phase space then yields 



NLO,l 
19 



v'^n'^ 1152 



^2 - 2 (1 - 



z In z 



2 In (1 - ^) + In ( ^ 



(164) 



Result 



In this subsection, we again summarise the results we obtain from the different sub- 
processes of gluon-induced Higgs production; we find total agreement with the results 
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in [591160]. We have 



_NLO,2 



a" 



486 V TT^v^ 



NLO,l I _NLO,2 



at 



576 V vr^f 2 



il-z) 



7-3z 



1 + ln 



^NLO,l I ^NLO,2 I ^ I ^ 



1 



576 V TT^f 2 

'ln(l - z 



6{l-z] 



11 



TT 



11 



+ 12 



1 - z 



2 V 6 3 

2 (2-2(l-;z)) ln(l-;z) 



2 H 



2 ^ ^ 



4.4 Higgs decay to gluons 

This process is used to validate the final-state collinear splitting functions 99,9 
q q, Eqs. fl82|) . fl7i|) . fl83|) . f l75|) as well as the final-final state interference term Eqs. flll4p . 
f lll5p . Besides the virtual correction to the vertex, both H ^ g g g and H ^ g q q 
contribute to the real radiation terms for this process. 



Leading-order contribution 

The leading-order squared and averaged matrix element for the process H 
d = 4 — 2e dimensions is given by 

integration over the two-particle phase space then yields 

r3 ^,2 



9 9 



e . 



(165) 



Virtual correction, charge renormalisation, and effective theory correction 
term 



The virtual matrix element is given by Eq. fll54p 

= I-Mlo 



2-K 



2 ( 4'Kfi'\ , / 1 , 2^2 



M; 



tt' + 0{e) X 2 
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which gives 

TvM = Flo ^ Ca ' r(l + e) (^-^ + ^vr^ + 0{e) ] . (166) 

The charge renormahsation term is given as in Eq. fll56p 

Teh = (4Z,)rLo(i^^ 

= -Flo 



r(l - £) \ fip J \ TT / e \ IT / /i 



/4 

,2 ■ 



(167) 



where we used that 1(1 - e)T{l + e) = 1 + 0{e^). We also need to add the NLO 
correction to the matching coefficient 

Teff = Tlo y . (168) 

Real emission 

For the real-emission contributions, both 

H ggg, H gqq 

need to be considered. The squared matrix element for H{Q) — > g{pi)q{p2)<i{Pi) is 
given by 



\M{H-^gqq)\ 



,2 16 al {pi + ps)^ + iPi + Ps)^ 



9 7rf2 (p2 +p3y 

16 al i^i + X2- If + (1 - X2f 
9 7rt>2 (1 — xi) 



where we used Xi = ^^^^ and with = = Mjj. The corresponding subtraction 
term is given by Eq. (174|) 

\ o (1 - 3^2) (1 - 0^3) 



1 « ^3 1 



9 7ri;2 \ — Xi 



x\ 



where the factor 2 arises because each of the gluons in the leading order contribution 
can split into a quark-antiquark pair. Note again that, in the case of non-unit jet 
functions, this factor needs to be replaced by 

2 ^ F(2)(p,J + F(2)(p,J (169) 

where F^'^\pi) denotes 

y 

Pi = Pi = P2 + P3 - Pi 

in the respective mappings. Integrating over the three-particle phase space and sum- 
ming over final state quark flavours yields 

dT{H ^ gqq) = ^ / ^^^3 (|-M(// ^ gqq)f - D,,,) = r^oy ("^^z) • 
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For the purely gluonic decay, we obtain for H{Q) — )• g{pi)g{p2)g{p. 



3J, 

^\MiH^ggg)f = 

32 1 r2s?23-5l2 , 2s?23Si3 , 2s?23'523 , 2si2Si3 2si2S23 , 2S13S23 , ^ 

H 1 1 1 h 0S123 



TTV^ 3 3! [ S13S23 S12S23 -312513 S23 Si3 S12 

where 



Stj = (Pt+PjY, S123 = {Pl+P2+P3f, 
S12 Si3 S23 

x = , y= , z = . 

-3123 ■3123 -3123 

In the subtraction terms, we have to sum over pj,j = 1, 2, 3: 

{Pj,Pi) = {PuP2), {Pj,Pe) = {pi,P3), {Pj,Pe) = {P2,P3), 
where we have the transformations 

. ^ . y . 1 . 

Pe = Pe+ Pj - Pfc, Pk = Pk 

1 1-1/ 

and we have to account for the fact that pk = {Pgi,Pg2) for each of the above scenarios. 
For the soft subtraction term, we therefore obtain 

n - ^ \^Tl^f\(^ ^^llA/f |2 - 32 1 [ 2g^23-3l2 2gj23-3l3 2g^23'323 

6 TIV 6 6\ I S13S23 S12S23 S12S13 

(170) 

where is given by Eq. flll4p and the sum runs over all possible combinations. 
The factor 1/3 in the above sum reflects the fact that the final-state gluons are indis- 
tinguishable!^, and we implicitly included a factor 2 in this expression which, in case of 
non-unit jet functions, needs to be modified according to Eq. fll69p . For the coUinear 



^In principle, the real emission matrix element includes all possible 3! combinations of {pi^pj^pk); 
for each of these, the matrix factorises according to 

Mrea.1 ^ V MlO 

in the singular limits. We obtain a symmetry factor ^ in front of the squared matrix element; however, 
the factorised splitting function is the same for pi ^ Pk, which effectively leads to a factor fr = -j- 
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contribution, we have 

= I J2'^999liPe^Pj)] \MlO 
4Si2Si3 



Snv"^ 3 1 S23(si2 + sn)^ 



[Sl2 + S23)(-Sl3 + -523) 



+ 



2s23g' 



{Si2 + S23)(Sl3 + S23) 



4Si2S23 



Sl3(Sl2 + 523)^ 
4S13S23 



2 - 



4Si3Q' 



Sl2(Sl3 + 523)^ 



,Sl2 + Si3)(Si3 + S23, 

4si2Q2 



+ 



+ 



S12 + S23)(Sl2 + S13) (Sl2 + S23)(Sl2 + S13 




where T^ggg is given by Eq. 

Integrating over the three-particle phase space yields 



■.NLO,3 
999 



(171) 



Integrated subtraction terms 



The two-particle phase-space contribution from the integrated gqq splitting function, 
Eq. (175]), yields 



^NLO,2 
999 



Flo -Tn 



1 



/ 47r/i 

7r^"r(l-.) [-Ml 



2\ e 



2 16 



where we have summed over final-state quark flavours and used 2pi ■ Q = Mfj for 
/ = 1, 2 being the gluon indices in the two-particle phase space. For ggg splitting, we 
use Eqs. (l83l) and (11151) for the integrated splitting function contributions to obtain 



■^NLO,2 
999 



27r V{l-e) V M- 



^ y 35 9 4 



(172) 



Result 

Combining all contributions to the total decay width at NLO, we have 
Tlo 
= Tlo 



■pvirt I ■pNLO,2 , ■pNLO,2 , ■pNLO,3 , -pNLCS 
~^ 9QQ ~^ 999 ^ giQ ~^ 999 



1 + — 

vr 



95 7 



+ Ech + Eeff 
2 



4 6 

which reproduces the result in the literature [63] . 



InfTA In^ 



H. 



:i73) 
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4.5 Deep inelastic scattering 



For a check of the (integrated) initial-final and final-initial interference terms given by 
Eqs. (11091) . (11141) . (IllUj) . and (I116p . we consider a subprocess of deep- inelastic scatter- 
ing: 

e~{Pi) + q{pi) — y e'{Po)q{p2) (174) 

at NLO, with a pure photon exchange in the t -channel. This suffices to test the validity 
of initial-final and final-initial state interference terms, therefore completing our check 
of all integrated subtraction terms for LO processes with up to two particles in the 
final state. 

4.5.1 Tree-level contribution 
Matrix element 

For the process fll74p . we have 

^El-^li°rn = 2e^QS^^, (175) 

where s = {pi+pi^ = {Po + P2)^, t = {Pi-PoY = {pi-P2f-, u = -s-t. Note that 
the singularity t = needs to be regularised by a proper definition of the jet function. 

Phase space 

The two particle phase space is given by 



4.5.2 Virtual matrix element 

The loop correction to the vertex is 

- %r r(lb) O' (-1 - ? - « ^ °(^) I ^-»- 

and we obtain 



\Mvr = 2Re{MBornM: 



virt) 



Svr r(l - e) \2pi-poJ [ e 
- |MBor„r|^ (-4 --(3 + 2 1nir-27)-8 - 3 (lnir-7) 



- W K + \n K + Y\ (177) 

6 



where K = 2tc fi"^/ {pi ■ po) 
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4.5.3 Real emission 



For a check of the remaining interference terms, the only relevant real-emission process 
we need to take into account is the quark-induced channel 

e'{Pi) + q{pi) — > (i'{Po)q{Pi)9{p-i)- (178) 

A complete discussion would also include gluon-induced processes; however, the cor- 
responding (integrated) subtraction terms Eqs. fl95l) . fl96|) have already been discussed 
in the context of single-W production in Section I4.1[ and are therefore not considered 
here. 



Matrix element 

The squared and spin-averaged matrix element for the real emission process is 



8 TT cts Ql 



spins 



{Pi ■ Ps) {P3 ■ Pa) {Pi ■ Po 

+ {Po ■ Plf + {Po ■ Pif] ■ 



Cf [{Pi ■ Pif + {Pi ■ Pif 



Three-particle phase space 

In our code, we follow the three-particle phase-space parametrisation of [51]. We then 
obtain for the outgoing four-vector^ 



P3 



Pi 



iPsI 



, Po 



where 



/ 1 \ 

sin^ 



\ cos6' / 

■5 - IpsI - Ipo 

P3 - Po 



COS^ 



IPol 



/ 1 \ 

COS 9 cos 7] sin ^ + sin 9 cos ^ 
sin 7] sin 

\ cos 9 cos ,^ — sin 6' cos 77 sin ^ J 



|P4 


2 _ 


|P3 


2 _ 


Ipo 


2 


2 1 


PsI 


IpoI 





P4 



Q-P3- Po- 



The three-particle phase space is 







dpi I dpi 



1 1-2-K 

d cos 9 I dip 
-1 Jo JO 



2tt 



(179) 



dr]. (180) 



^We are aware that a more detailed analysis could be optimised by a more process-specific 
parametrisation of phase space which specifically maps the singularity structures, in combination with 
a multi-channel integration as e.g. VAMP [65'; however, for our purposes this simple parametrisation 
proved sufficient. 
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4.5.4 Subtraction terms 
Initial-state subtraction 

The initial-state subtraction is given by Eqs. (1^ and fll09p : 

= Ci. 1 - a; - y + — ^ -MBorn P ' 

xypi-Pi V v[l-x) + yj 

where lA^BomP is given by Eq. fll75p and with 

Po-pA Pl-Pz ~ Pi -Pi {Pl ■ Pi) {P^ ■ Pi) ..g.x 

X — y — z = V — (lolj 

Pi-Pi Pi-Pi' (P4 ■ Q) (P3 ■ Q) ' 

The final-state momenta p4, Po are mapped according to the corresponding transfor- 
mation specified below. 

Initial-state kinematics: momentum mapping 

The momentum mapping needed here has been specified in Section l3.1.2j for the initial- 
state momenta, we obtain 

Pi = Xpi, Pi = Pi. 
The four- vectors of the outgoing particles are transformed using {w = o, 4) 

= A^,{K,k)pl, 
where A(i^', K) is defined according to Eq. f lT7|) and with 

K = xpi +pi, K = pi+pi-ps. 

Final-state subtraction 

The final-state subtraction term is given by Eqs. fl7T]) and flll4p 



D 



4 3 _ AnagCp 
y {pi ■ Pi) 

where we use 



y J, ^ {l-v){l-z{l-y)) 



1-y V [l-z{l-y)] + y[{l-y)a + l] 



|-MBorn(p)P 



A . . PS-Pa - Po 

P4-Q = Pi-Pi, y = - — ^ 



p3-po p3-Q ~ 1 Pl-Po 



X 



p3-po+Pi-po' p3-Q+p4-Q' 1 Pi -Ps +P1 ■P4' 

Pi ■ P3 „ ^ {p3 ■ Po) {Pa ■ Po) 



P1-P3+P1-P4 {p3-Q)'^ 

The momenta in the Born matrix element need to be calculated using the transforma- 
tion defined below. 
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Final-state kinematics: momentum mapping 

By definition, we have 

1 



Pi = Pi, Pi = Pi, Pa 



[Pz + Pa - y {Pi + Pi)] ■ 



Po 



Energy and momentum conservation then implies 

Po = Pi+Pi- Pa 

4.5.5 Integrated subtraction terms 
Initial-state subtraction 

The initial-state integrated subtraction term, including the interference term, is given 
by Eqs. f l93p and flllOp respectively, and we obtain 



Pie) 



KHx:e) 



1 1 /3 , . 



1 

X 



e \2 
1 fl + x' 



i Li2 (1 - 5o) + ^ In' 5o + 2 ll,{~Zo), 



e \ 1 — X 

'ln(l — x) 



'l-x)2 Ind -x) 



l+x^ 
1 — X 



In X 



+ 4x 



1 — X 



:i82i 



where z is defined in Eq. fllSip . The four- vectors in z now have to be calculated from 
{p} using the inverse mapping 



Pi = -Pi, Pi = Pi, Pa = A{K,K)pi 

X 

where A is defined in Eq. (|T7|) and 

X + y'{l — x) 



[183) 



K = pi + Pi, K 



X 



-pi + [l-y'{l- x)] Pi + \k_i\k_i 



in terms of the m-particle phase space variables. Here, 



:i-x) 



:i - y') y' 

X 



2piPi, k_i 



( \ 

2f - 1 



v 



-2jv(\-v) 




The complete subtraction term is given by 
1 /47^/i2^" 



C 



2tx ^ r(l -e) \2pi-p 



dx 



5{x-l)r{e) + k\x]e) + P^o\\{x,e]^F) lA^Bom 



(184) 
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with 



Pcon{x,e; jj^p) 



11 f 2pa -PbY A + a; 



I \ \ — X 

The integrated subtraction term contains two finite integrals 



C(5o) 



dy 



y' v/4y'2(i_5o) + 52 
X In 



2v/4y'2(l-5o) + 5o'\/(W) 



i r ^ 

(l-x)+ ttJo 



2 1/' + zo - 2 5o + v/4y'2 (1 - zo) + A 

- 1 



v^t; (1 - V) N{x, y', z, v) 



(185) 



which need to be integrated numerically. Here, 

Pa-pi 1 



N 



pi-Q^-x 



y 



and p3 needs to be reconstructed using Eqs. (185|) . (186|) . and (l87|) . For the implementa- 
tion of the + distribution, we also have to use 



No = N{x = 1) = -T, h y . 

Pi ■ Q 

If we expand in e, we then obtain for the subtraction term 

a 



27r 



dx < 5(1 — x) 



h^\v2~ ^""^^ + In A - 7j + ^Li2 (1 -5o) 



+ ^ In'io + ^In'A - 7ln A-^ + + ("^ - In 5o ) (In A - 7) + 2/fi„(5o) 



1 

+- 

a; 



2 (1 - x) In (1 - X 

1 + x^" 



+ 



1 — X 



In B 



l + x^\ , , /In(l-x) 

lnx + 4x ^ - 

1 — X I , \1 — X 



+ iL{~z,x) \\MboAp)\' 



with 



A 



27r/x2 
Pi - Pi' 



B 



2pi • Pi 

1 ' 



and 7 being the Euler number. 
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Final-state subtraction 

The final-state integrated subtraction term is given by Eqs. 0721] and 01161] : 



1 1 

- + - 

e 



+ \ia{do + 1) 



1 + y - 2 In (2) ln(ao + 1) + 2 l^'(«o + 1) 



ao 



where 



do + 1 

1 du 
u 



1.. 

2 ^ 



2 In 2 + 



ao 



(3-0 + 1 



+ C(5o) + IIM 



X In 



du dx 



V^l + 4 5o(l + 00)^2 
(1-n) 



X 



(^l + 2 5oU + v^l + 4 5o(l + So)u2y 

x(l— x + MX [(1— Mx)a + 2]) 
k{u, X, a) 

1 



^l + 4ao M2(i + ao) 



(186) 



Here, 



k^{x, u, a) = [(1 + ux — x){z — z') + ux ((1 — ux) a + 1)]^ 
+ 4:ux z {1 — z) (1 + ux — x) ((1— 'ux)a + l) 



and 



xil — u) , _ _ 

z = , z = uxa, a 



Pi -Po 



1-ux ' ' {l-y)pi-Pi + ypi-Pi 

in terms of the m-particle phase-space variables. The complete subtraction term is 
given by 

■27r/i2\^ 



a 



1 



2 7r r(l -e) \pi-p 
Expanding the complete integrand in e, we obtain 

a 



lf{s)\MBoUp)?. 



2tx 



e 

1 - .2,. 



ln(ao + 1) + In D - 7 



TT 



1 + — - 2 In (2) ln(So + 1) 



+ - ln"(ao + 1) + -Lis 



ao + 1 



Lio 



ao 



0,0 + 1 



In^ L> + 72 



7 In D 



+ - + ln(5o + 1) (In - 7) + C(ao) + ll^{d) )\MBorn 



with D = 2 TT fi"^ / {pi ■ pi) . 
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Effective /, K, P terms 



We can combine tlie initial- and final-state integrated subtraction terms to an effective 
/, K, P term which should be applied in the spirit of an initial-state subtraction term 
(i.e. the / term is multiplied with 5 (1 — x)). We then obtain 



a 



I i.^) = TT^Cf X <^ ^ + - [3 - 2 In^o + 2 In A - 27] 



27r 



1 



+ 3 Li2 (1 - zo) + 2 In 2 In zo - 2 [(1 - ^o)'] + In^ zo 



+ (In A -7)' 



TT 



(3 - 2 ln5o) (In A - 7) - 1 + C '°(^o) + iLC^) 



a« ^ I 1 
27r \x 



2(l-a;) In(l-x) - 
ln(l — X 



1 + x^ 
1 — X 



In X 



a. 



+ 4x 
1 fl + x^ 



1 — X 



2 vr X \ 1 — X 



In 



2pi • Pi 



(187) 



where 



Y tot,0/ ~ 
'fin 



(5o) 



' dy 



y 1 v/4y2(i-5o) + 52 



X In 



2z - y) 



In 2 



(2y + Zo-2yzo + ^Ay^ (1 - Zq) + ll ^ 

/gj^(a) is given by Eq. fll86p and /g^(^,x) by Eq. (11851) . We can further simplify 

'27r^2^ 



In ( - 

,^0 



In 



Vi-Vo 



Combined two-particle phase-space contribution 

Adding the /, P terms from the last section to the virtual contribution, we have 



dx\M\ 



dx 



as 
27 



Cf5{1-x) 



-9 



- \u,[{l-~zof] 



+ 2 In 2 ln5o + 3 Inio + 3Li2(l - 5o) + C'°(^o) + llM 
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4xl0'5 



3xl0'^ 



2xl0'"^ 



IxlO'^ 




Figure 6: NLO (green, dashed) and LO (blue, solid) partonic cross sections for DIS 
subprocess eg — )■ eq{g), as a function of parton level (HERA-like) cm energies, with 
angular cuts cos 6ee < 0.8. Partonic cross section prior to convolution with PDFs. 
Relative NLO corrections are around 3.4% 



where K^^{x;z) and Pg°*(x;/i|') are given by Eq. f ll87p . Note that ll^{z,x) contains 
the four- vectors p4 and ps, which need to be reconstructed using Eqs. f ll83p and (185|) . 
dSSD, dHZD, respectively. 



4.5.6 Results 

We have compared the above results numerically with an implementation of the Catani- 
Seymour dipole subtraction; the corresponding terms can easily be obtained from [31] 
and are therefore omitted here. Figs. Eland [7] show the behaviour of the total partonic 
cross sections and the differences between the application of the two schemes at par- 
ton level for varying (HERA-like) center-of-mass energies, where we applied a cut of 
cos 9ee < 0.8 to cut out terms for which t = {pi — PoY = 0. We see that the results 
agree on the per-mil levelll, therefore verifying the nontrivial check of our subtraction 
prescription with a modified final-state mapping. In Fig. |H1 we show the behaviour of 
the differences between the two- and three-particle phase-space contributions for vary- 
ing center-of-mass energies, again verifying that their cancellation is non-trivial as the 
contributions from different phase space integrations vary widely in magnitude for the 
two schemes. The results for our new scheme have been obtained using subtractions in 
the vicinities of the singular regions in phase space onlj|§. A more detailed investigation 
for this process, including the implementation in a parton level Monte Carlo generator, 
is in the line of future work. 



^Note that this result has been obtained with a relatively mild angular cut. 

^This idea has been documented in [3D] in the context of FKS subtraction and in \EW for Catani 
Seymour dipoles. Explicit expressions for our scheme will be presented elsewhere. 
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Figure 7: As Figure El relative difference between NLO contributions using Nagy-Soper 
(NS) and Catani Seymour (CS) subtraction terms. Errors are integration errors; results 
agree at the permil-level. 
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Figure 8: As Figure [61 Behaviour of the difference A ' 



NS,CS 



cr 



(2,3) 
NS 



,(2,3) 



— cr^q ' between 



the two schemes for the two particle final state phase space (green, dashed) and three 
particle final state phase space (blue, solid), respectively. In the sum A^g + A*-^^ 



NS,CS 



(black dots), the large differences cancel. In the new scheme, subtractions have been 
restricted to singular regions. 
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5 Conclusion and Outlook 



In this paper, we propose a new NLO subtraction scheme which is based on the mo- 
mentum mapping and on the sphtting functions derived in the context of an improved 
parton shower formulation |l0]. One essential feature of our scheme is that we use a 
global momentum mapping in which all of the partons participate. As a result, the 
number of momentum mappings needed to evaluate the subtraction terms and the 
corresponding leading-order matrix elements is significantly smaller than for standard 
subtraction schemes. A further important feature of our scheme is that we derive the 
subtraction terms from the splitting functions that describe an improved parton shower 
with quantum interference. The use of the shower splitting functions as subtraction 
terms greatly simplifies the matching of a NLO calculation with the corresponding 
parton shower. 

To establish our scheme we have focused on the simple case of collider processes 
with up to two massless particles in the final state. We have presented formulae for 
all subtraction terms and integrated counterparts needed to address such processes at 
NLO, and have applied the results to a variety of basic lepton- and hadron-colhder 
processes. In all cases, we have reproduced the results from the literature and have 
shown that our implementation agrees with results obtained using the Catani-Seymour 
dipole subtraction. In the scheme proposed here, the mapping of parton momenta is 
given by a generic description for initial- or final-state emitters. The finite contributions 
in the integrated soft interference terms between initial- and final-state partons thus 
depend on the kinematics of the process and need to be integrated numerically. The 
implementation of initial-final and final-initial interference terms therefore accounts for 
a non-trivial check of our subtraction prescription. 

The advantages of two main features of our scheme, i.e. the global mapping and 
the use of subtraction terms derived from the splitting functions of an improved parton 
shower, become most apparent when applying the scheme to NLO multi-parton pro- 
cesses and/or matching the NLO calculation with the corresponding parton shower. 
Both applications are in the line of future workjf] In this paper, we have demonstrated 
in a first step that the splitting functions of the improved shower, in combination 
with the corresponding global mapping, can be used as local counterterms in a sub- 
traction scheme for processes with relatively simple final-state kinematics. While the 
subtraction terms for the real-emission contributions can readily be applied to generic 
multi-parton final states, the mapping prescription for final-state emitters leads to 
more involved finite parts of the integrated subtraction terms when considering pro- 
cesses with three or more final-state particles. Although the integrated subtraction 
terms for the general case have been privately available for some time, and an example 
for the final-state splitting function g ^ qq has been presented in [H], a generic appli- 
cation to a more challenging physical process is still work in progress. However, we are 
confident that results for multi-parton final states, which will allow for a more powerful 
test of the features of the new scheme, will become available in the near future. 

^We note that work is underway [67] to implement the scheme presented here into the Helac Event 
Generator framework [11) . 
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£ fi fi fj Vi X J: — colour 
F q q g £f,{Pj,Sj]Q) ' . ^ 



2pe-ni [iPe-Pj] 



I q q g -e^[pj,Sf,Q) ^ ^ ^ 



1 g q q —e'^[Pi,sg,(^) Ufj,u[Pe — Pj,ne) — — ^ t 



U {pj, Sj)-1>'[f^ - p^]T^iV{ pi, Si) 

2pi-ni [{pe-pjy^ 



I q g q -ef,{pe,Si;Q) tr- - ^^?^ ^ 



Table 2: Splitting amplitudes vi{{p, f}j, Sj, si, si) involving a qqg splitting. We have 
removed a common factor a/IttoJ removed. The label i denotes either initial-state 
indices / = {a, b} or final-state indices F = {l,...,m}. The light-like vector ne is 
defined in Eq. fllQip . Taken from |40] . 

A Splitting amplitudes 

In Tabled we list the splitting amplitudes for qqg splittings as given in [10]. For triple 
gluon splittings, we have for the final state 

M{P, f}m+l,Sj,Se,Si) 



EaiPj, Sj] Qyep{pi, Si; QYe^ipe, se; Q) (188) 

Zpj-Pi 

X v°''^'^{pj,pi, -pj - Pi) D^y{pi + pf, ne) . 



For an initial state splitting, we have 

ViiiP, f}m+l,Sj,Si,Se) 



^'^^^ £a{Pj, Sj] Q)*ei3{pe, Se; Q)e''{pe, se; Q)* (189) 



2pj -pe 

X v"''^^ (pj , -pe, Pe - Pj ) D^^ {pe - pj ; ne) . 

We use standard notation where U{p, s), U{p, s), V{p, s), V{p, s) denote spinors of the 
fermions with a four-momentum p and spin s, and £q(p, s; Q) are the gluon polarisation 
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vectors. The ggg vertex has the form 

v"''\Pa,Pb,Pc) = g"^{Pa - PbV + g^\pi, - PcT + g^'^iPc - Paf . (190) 

The transverse projection tensor D^^^p^ ~Pj'^ ^e) is defined according to Eq. fHTl) . The 
hght-hke vector rii is given by 



Pb , £ = a , 

Pa , i = h , 

Q-p^ + viQ-Per 



(191) 



B Integration measures in terms of singular vari- 
ables 

In this section, we give the integration measures for the initial- and final-state splittings 
in terms of the singular variables, and relate these to the four-momenta in the real- 
emission phase spaces. 

B.l Initial-state integration measure 

The initial-state integration measure in d = 4 — 2e dimensions is given by Eq. 
and related to the variables in Section 13.3.21 via 

= dxdy' f(^°4)^(7vr)^'-' ~ "^^'"'^ " ^'^^^^ ® " ® " ^'^ ^'^ ' 

(192) 

where the center of mass energy in the m + 1 phase space is given by s = rjaVb^ = 
{2pa ■ph)/x. Here, a; — >• 1 corresponds to the soft and — )■ to the collinear singular 
limit. For the integration of the interference terms, we additionally have to parametrise 
the azimuthal angle of pj in the integration measure; this parametrisation is frame 
dependent. We obtain 



di, = dxdy'dv ' ^^r^^ — ^rp^(i-^)'"'Mi/'(i-i/'rM^(i-^r^ 

X e [(1 - x) x] [(1 - y') y'] e[v{l- v)] (193) 

where v = | (1 — cosy?) is related to the azimuthal angle of pj in the center for mass 
frame of Pa,Pb and pk defines the x, z plane. 

B.2 Final-state integration measure 

The final-state integration measure in ci = A — 2e dimensions is given by Eq. fl68|l . In 
terms of integration variables which parametrise the singularities of the integrands, we 
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obtain 



x(2 7r)'^5(^) (P, + P^-(l-y)p,-yQ) 
i2p,-Qy-^ (Any ri 



I duu-^il-u)-' f -a;)-^(194) 



16 7r2 T{l-e) 

where u and x are related to the m + 1 particle kinematics via 

Pj -Q Pi- Pi 

Pe-Q Pj -Q 

Here, x — )■ corresponds to the soft and m — )■ to the collinear limit of the integration. 
While the above parametrisation suffices in the integration of the collinear subtraction 
terms, we need to introduce an additional angle between emitted parton and spectator 
in the interference terms. In contrast to the above parametrisation, the additional 
integration variable is frame dependent; in the center of mass system of pe,pj, where 
pk defines the x, z plane, we can write 



16 r(i-.) 



1 r'X A-l 



X / duu-^{\-u)-^ / dxx^-^^ {\- x)-^ / dv\v{\-v)\ ^,(195) 
Jo io io 

with V = I (1 — cos(f) is related to the azimuthal angle of pj. The above parametri- 
sations were used in all integrations presented in this paper. 



C Altarelli-Parisi splitting functions 

This section contains a list of the well known Altarelli-Parisi splitting functions |15] . 
which are evolution kernels of the DGLAP equation [691 [70lll5l [71] . They also describe 
the behaviour of parton splittings by giving the probability of finding a parton of type 
b with momentum fraction a; in a parton of type a in the collinear limit: 



a{p) b{xp + k^+ 0{kl)) + c {{l-z)p-k^+ 0{kl)) 
At leading order, the splitting functions are given by 



(196) 



P 



X] 



l + x2 3^,^ 



X] 



P33(x) 



Tr [x^ + (1 - Xf] , Tr 



Cf 
2Ca 



l + (l-x) 

X 
X 



+ 



X 



XI 



X 



+ X {1 — x) 



+ 6(1 



X 



llCA-AufTn 
6 



(197) 
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where Uf is the number of quark flavours in the theory. The + distribution is defined 
in the standard way 

1 pi pi pi 

f{x)g+{x)dx= / g{x){f{x) - f{l))dx = / g{x) f{x) dx - f{l) / g{x)dx 
Jo Jo Jo 

(198) 

for the convolution with a test function f{x). 
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